
Calhoun 

Iniiiiuiiortfl Arthivcof (he Navjl Pwigndualt School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2013-12 

Automated techniques for rapid analysis of 
momentum exchange devices 

Smythe, Reid W. 

Monterey, California: Naval Postgraduate School 


http://hdl.handle.net/10945/39014 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


htt p ://w w w. n ps.e-du/l ib ra ry 


Calhoun is the Naval Postgraduate School's public access digitaI repository for 
research materials and institutional publications created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed —and published —scholarly author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 







NAVAL 

POSTGRADUATE 

SCHOOL 

MONTEREY, CALIFORNIA 


THESIS 


AUTOMATED TECHNIQUES FOR RAPID ANALYSIS OF 
MOMENTUM EXCHANGE DEVICES 

by 

Reid W. Smythe 
December 2013 

Co-Advisors: Mark Karpenko 

Isaac M. Ross 


Approved for public release; distribution is unlimited 




THIS PAGE INTENTIONALLY LEFT BLANK 



REPORT DOCUMENTATION PAGE 


Form Approved OMB No. 0704- 
0188 


Public reporting burden for this collection of information is estimated to average 1 hour per 
response, including the time for reviewing instruction, searching existing data sources, gathering 
and maintaining the data needed, and completing and reviewing the collection of information. Send 
comments regarding this burden estimate or any other aspect of this collection of information, 
including suggestions for reducing this burden, to Washington headquarters Services, Directorate 
for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 
22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188) 
Washington DC 20503. 

I. AGENCY USE ONLY (Leave blank) 2 . REPORT DATE 3. REPORT TYPE AND DATES COVERED 

December 2013 Master's Thesis 

4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 

AUTOMATED TECHNIQUES FOR RAPID ANALYSIS OF MOMENTUM 
EXCHANGE DEVICES _ 

6. AUTHOR(S) Reid W. Smythe 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING ORGANIZATION 

Naval Postgraduate School REPORT NUMBER 

Monterey, CA 93943-5000 

9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING/MONITORING 
N/A AGENCY REPORT NUMBER 

II. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and 
do not reflect the official policy or position of the Department of Defense or the U.S. 
Government. IRB Protocol number_N/A. 

12a. DISTRIBUTION / AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Approved for public release; distribution is unlimited 

13. ABSTRACT (maximum 200 words) 

This thesis develops a suite of automated techniques aimed at analyzing the 
characteristics of arrays of control moment gyroscopes (CMGs) and/or reaction 
wheels (RWs). Three specific areas, relevant to the design of new spacecraft, are 
examined: momentum space, torque space, and reaction forces. The momentum space 
analysis creates a maximum saturation envelope for both pure CMG/RW arrays and for 
hybrid arrays. The torque space analysis creates a maximum envelope for RWs and 
presents the idea of state-space analysis for CMGs, which leads to the ability to 
determine singularity free operational envelopes. These envelopes allow satellite 
slew performance to be estimated, allowing for an initial determination of 
component size. The reaction forces and moments seen at the attachment points of 
the RWs and CMGs during maneuvers is also analyzed. This analysis, which has not 
been reported in the literature, utilizes telemetry data and the geometry of the 
spacecraft to explore the relationship between maneuver types and the loads seen by 
the CMG/RW mechanism. 

14. SUBJECT TERMS Reaction Wheel, Control Moment Gyroscope, Reaction 15. NUMBER OF 

Forces PAGES 

_169_ 

16. PRICE CODE 

17. SECURITY 18. SECURITY 19. SECURITY 20. LIMITATION OF 

CLASSIFICATION OF CLASSIFICATION OF THIS CLASSIFICATION OF ABSTRACT 

REPORT PAGE ABSTRACT 

Unclassified Unclassified Unclassified UU 

NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Std. 239-18 


1 




























THIS PAGE INTENTIONALLY LEFT BLANK 


11 



Approved for public release; distribution is unlimited 


AUTOMATED TECHNIQUES FOR RAPID ANALYSIS OF MOMENTUM 

EXCHANGE DEVICES 


Reid W. Smythe 

Lieutenant, United States Navy 
B.S., United States Naval Academy, 2007 


Submitted in partial fulfillment of the 
requirements for the degree of 

MASTER OF SCIENCE IN ASTRONAUTICAL ENGINEERING 

AND 

ASTRONAUTICAL ENGINEER 


from the 


NAVAL POSTGRADUATE SCHOOL 
December 2013 


Author: Reid W. Smythe 


Approved by: Mark Karpenko, PhD 

Co-Advisor 


Isaac M. Ross, PhD 
Co-Advisor 


Knox T. Millsaps, PhD 

Chair, Department of Mechanical and 

Aerospace Engineering 

iii 



THIS PAGE INTENTIONALLY LEFT BLANK 


IV 



ABSTRACT 


This thesis develops a suite of automated techniques aimed 
at analyzing the characteristics of arrays of control 
moment gyroscopes (CMGs) and/or reaction wheels (RWs). 
Three specific areas, relevant to the design of new 
spacecraft, are examined: momentum space, torque space, and 
reaction forces. The momentum space analysis creates a 
maximum saturation envelope for both pure CMG/RW arrays and 
for hybrid arrays. The torque space analysis creates a 
maximum envelope for RWs and presents the idea of state- 
space analysis for CMGs, which leads to the ability to 
determine singularity free operational envelopes. These 
envelopes allow satellite slew performance to be estimated, 
allowing for an initial determination of component size. 
The reaction forces and moments seen at the attachment 
points of the RWs and CMGs during maneuvers is also 
analyzed. This analysis, which has not been reported in the 
literature, utilizes telemetry data and the geometry of the 
spacecraft to explore the relationship between maneuver 
types and the loads seen by the CMG/RW mechanism. 
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I. 


INTRODUCTION 


Changing the attitude of an object generally relies on 
a common technique: modifying its angular momentum. The 
mechanisms available to effect this change differ depending 
on whether the object resides in a medium such as air or 
water or if it exists in a vacuum. Everyday experiences 
such as the turning of a car or boat or the banking of an 
airplane have given many people an intuitive feel for how 
objects react on Earth. This intuition does not translate 
to objects in space because the motive forces used on Earth 
are not available: there is nothing in the environment to 
push off of, hold on to or generate lift from. One 
solution that has been developed is to trade the angular 
momentum of the spacecraft with a component mounted inside 
the spacecraft. This allows the total angular momentum of 
the system to nominally remain constant, conservation of 
momentum, while still allowing the spacecraft to spin via 
momentum transfer. These devices are generally referred to 
as momentum exchange devices (MEDs). 

A. MOMENTUM EXCHANGE DEVICES 

Although there are several different types of MEDs, 
the majority of those used today can be divided into two 
groups. The first group, reaction wheels (RWs) consist of 
a relatively large mass (a rotor) mounted on a shaft and 
attached to a motor. The motor allows the RW to change its 
angular momentum, by changing the rotational speed of the 
rotor, and thus the angular momentum of the spacecraft to 
which the RW is mounted. The second group of momentum 


exchange devices is single gimbal control moment gyroscopes 
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(CMGs) . The basic idea behind CMGs is to mount a reaction 
wheel rotating at a nominally constant speed on a second 
shaft that can gimbal. As the second shaft gimbals, the 
direction of the momentum vector of the rotating wheel 
changes, thus changing the angular momentum vector of the 
spacecraft as a whole. Note that variable speed CMGs can 
be implemented to combine the characteristics of both RWs 
and CMGs into one device [1, p. 179]. 

Combining several MEDs in different orientations 
results in an MED array and allows control of the 
spacecraft momentum along additional dimensions. 

B. OBJECTIVES AND SCOPE OF THESIS 

When conducting an initial design of an attitude 
control system (ACS), attitude control devices are selected 
based on the mission of the satellite, gross design 
characteristics (spin-stabilized, tumbling, three-axis 
stabilized) , mission life and other big picture 
specifications. Using back-of-the-envelope calculations, 
initial sizing is performed, and then a simulation is 
developed to verify and iterate on the decisions. This 
process continues until a final design has been established 
that meets various requirements such as mass, performance 
and service life. 

The goal of this thesis is to develop equations and 
algorithms to help automate the process of analyzing the 
expected performance of an arbitrary MED array. This will 
assist in determining an initial component size and 
location more accurately than back-of-the-envelope 
computations, based on required performance. 
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The developed algorithms are designed for ease of 
implementation in a vector- or matrix-based language such 
as MATLAB and can utilize the benefits of parallel 
processing when available. 

The initial analyses will consider the well-studied 
areas of angular momentum and torque. These provide 
important characteristics of the array performance. By 
using a clearly defined system to describe the orientation 
of MEDs inside a spacecraft, universal base equations were 
derived to eliminate the need for the creation of one-off 
solutions for a given array. These equations are not 
restricted by array size or component orientation. 

One issue with performing a direct comparison between 
RWs and CMGs comes from the fact that any torque trajectory 
for a system of RWs results in a continuous trajectory in 
state-space. This relationship stems from the fact that 
torque is not state dependent. This is not the case for 
CMGs, as the torque available to a CMG at any given moment 
depends on the current state of the CMG array. To overcome 
this problem, a state space analysis for CMG arrays is 
developed as part of this thesis. 

The final objective of this work is to examine the 
reaction forces and moments seen by the MEDs as a 
spacecraft maneuvers. This particular area has received 
little to no attention in the literature. 

In general, momentum exchange between MEDs and the 
host spacecraft is discussed in precisely that term: 
angular momentum of the system as a whole is conserved by 
exchanging angular momentum between the two bodies. The 
mechanism for this exchange, however, has a physical 

3 



consequence in the form of applied forces and moments on 
the attachment points between the spacecraft and the MED. 

When operating under rigid body assumptions and 
modeling spacecraft maneuvers, these reactions can be 
ignored as they are conservative. However, they represent 
possible failure modes and thus must be either considered 
when designing the connections, such as bearings, otherwise 
the connections may be significantly over-engineered adding 
to cost and mass. 

The force analysis will also allow for the comparison 
between different control laws. This can be useful, for 
example, if it is desired to extend the life of components 
by minimizing the peak forces and moments applied during 
maneuvers through a redesign of the attitude control 
system. 

C. THESIS LAYOUT 

Chapter II presents background material necessary to 
understand the derivations of equations in matrix form 
rather than vector form. 

Chapter III describes the derivation of the equations 
used to calculate the angular momentum space for RW and CMG 
arrays. It also presents the momentum inversion algorithm 
that can be used to find the gimbal angles necessary for a 
CMG array to orient the maximum net angular momentum vector 
in a given direction. 

Chapter IV derives the equations used to determine RW 
array torque space. For CMGs, it shows the background and 
principles behind state space analysis, including the 
determination of singularity-free momentum space. 
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Chapter V steps through the derivation of the 
equations used in the reaction force analysis. This 
includes determining both the inertial motion of the MED as 
well as the applied forces and moments used in the force 
and moment balance equation. 

Chapter VI presents several applications of the topics 
discussed in previous chapters. These applications include 
both notional array constructs and telemetry from 
previously validated simulations. 

Chapter VII gives the overall conclusions noted from 
the research and provides some areas that would benefit 
from future analysis. 
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II. MATRIX-BASED SOLUTIONS 


Although coordinate transformations and other vector 
mathematics associated with space systems engineering tends 
to favor vector or dyadic based solutions, the equivalent 
matrix form will be used for the derivations in this 
thesis. The primary purpose for this decision is to allow 
for simple implementation in matrix-based languages as 
MATLAB. The bulk of the following derivation is based 
heavily on Hughes [2, pp. 22-24]. 

Equation (1) gives one of the most used relationships 
when analyzing spacecraft dynamics: the transport theorem 
[1, pg. 12] . It allows a vector expressed in one frame to 
be differentiated with respect to another frame. 

—( a v) = —( a v) + co a,N x a v (1) 

Equation (1) is generally used to find the inertial 
derivatives of vectors that are more easily expressed in a 
non-inertial frame. The equivalent matrix expression. 
Equation (2), involves transforming the vector into the 
appropriate frame, then using the product rule to 
differentiate. 


N v= N a R a v 

N v= N R a v+ N R a v 


( 2 ) 


^R is the direction cosine matrix (DCM) that transforms the 
vector a v from the a frame to the N frame. 


Although Equations (1) and (2) are mathematically 
equivalent. Equation (2) is not written in a generally 
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useful form as the time derivative of the DCM between two 


frames is not a physically measured quantity. This results 
in additional computational overhead as that quantity must 
be calculated. To overcome this, the DCM rate can be 
rewritten as the product between the DCM and a skew- 
symmetric matrix of the rotation rates between the frames. 

First, consider the non-collinear frames a and N 
fixed in inertial space and described by the vectrices 1 [2, 
p. 8] {«} and {£>} , defined as Equation (3) . 



r /v h 

a x 


r ^ h 

II 

jaP 

a 2 

ll 

b 2 


a x 

< J J 
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The DCM ( R) between the two frames has the properties 
noted in Equation (4) . 

M = ;R{b] 

w-:*w 

lR = IR t = {/>} [a) T 

b n a n _ a n b n _ j 

a^b^ — ~ *3x3 

Now, allow the frames to rotate at different rates 
with respect to the inertial frame. These rates can be 
represented by the free vectors and . The 

relative rate between the frames is simply the vector 
difference between the two as given in Equation (5). 


1 The term vectrix (plural: vectrices) was "coined" by Hughes to 
describe a convenient notation with which to keep track of the 
multitude of reference frames used in spacecraft attitude dynamics. 
They are particularly useful when converting vector equations to their 
scalar equivalents [1, pg. 522]. 



alb a/Inertial b/inertial 

CO = CO —co 

,-bla 

CO = —co 


Equation (6) differentiates {a} in frame b with 
respect to time using the transport theorem and gives the 
expected result. 


v dt j 


d {a} 'j f d {c/} 
v dt 


+ co alb x l a 


a ( 


d {c/} 
v dt j 


= 0 


Differentiating {£>]■ in its own frame with respect to 
time yields an important equality, shown in Equation (7): 


Y 

V 


d { b } 
dt 


bf d h R\ , b h (dW' 

W\+ a R 


J 

0 

0 


dt 


V 


dt 


J 


= :*h + :*K«h) 

=(*«+.**<){«}, 


:7) 


= 


0 




alb 


-CO - 

0 




CO, 


-co. 


alb r^a/b 

CD, 


-CD X 

0 


alb 


Since | a] cannot be zero, this forces the parenthetical 
sum to add to zero to satisfy the equality in Equation (8). 


0 = h R + b Rco 

a a 


alb 


co a ' N = - b R b R T = b R b R T 


Returning to v, consider v as a free vector in R 3 
(Equation (9) ) . 


v = 
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To express the vector in either frame, use the inner 
product between the vectrix and the vector, given in 
Equation (10). 


v ={a} r a v 
v ={ b} T h v 


( 10 ) 


This changes Equation (2) to the form found in 
Equation (11). 


b 


d a v 

dt 


d{a] T a 


dt 


v + {a] T a v 


- +a"' , x{a} % 


dt 

dv a iu a 
— + co X a v 
dt 


( 11 ) 


Equation (11) matches Equation (1), which is an 
expected result. The penultimate step of Equation (12) 
involves expressing co alN in terms of frame a and then using 
matrix identities to rewrite Equation (11) in the same form 
as Equation (2) . 


a/N ( 1 T a a/N f a a/N~\ T ( ) 

co ={a| co = co {a} 


r d 

{a} 

T 1 

^a/N 

( \T } 

r a a/N ( 1 1 

V 

dx 

2 

= CO 

x{a} =( 

L ® J l a }} 


M t a a/N ( ) ( 1 T ( ) T a a/N 

co {a}(a| ={a\ co 


M x =M 


( 12 ) 


Finally, combining Equation (10)-(12) gives the 
desired result in Equation (13) . 

{b} T b v = {a} T a v + {a} T [“co alb Yv 
fc v = {iV}{a} r a v + {b}{a} T [ a co a,b ] a v 
b v= b R a v+ b R[ a <p alb ] a v 
10 


(13) 



By inspection with Equation (2), it is clear that 
Equation (14) holds. 

N a R = N a R a co a,N (14) 

A more detailed proof of this fact can be found in 
Likins [3, pp. 107-110] . 

Although performing second order derivatives on 
equations with several rotations can become unwieldy based 
on the number of terms generated with the product rule, the 
relatively small number of coordinate frames to be analyzed 
circumvents this problem. An additional benefit is seen in 
the removal of the cross-product as a function. This 
eliminates the need for additional function calls within an 
equation, improving both readability and speed. It will be 
necessary, however, to create the skew-symmetric matrices 
manually, slightly adding to the length of the code. 
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III. ANGULAR MOMENTUM SPACE 


The momentum space available to an MED array provides 
an indication of how much momentum the array can exchange 
with the spacecraft in any given direction, stemming from 
the fact that by the law of conservation of angular 
momentum. Equation (15), angular momentum must be conserved 
in the absence of external torques. 

h to tai = h s/c + h MED = const. (15) 

Physically, this "maximum momentum exchange" value 
represents the maximum rate the array can cause the 
spacecraft to spin (or the maximum spacecraft rate that the 
array can stabilize), shown by Equation (15). The shape of 
the MED array momentum space is unique to each individual 
array and is based on both the orientation and the physical 
and mechanical properties of each MED. 

h + h = I m +h MR) 

"'S/C "MED J SIC^S/C "MED \ ± ^ ) 


A. REACTION WHEELS 
1. Overview 

Consider a perfectly constructed RW rotating about its 
z-axis, as shown by Figure 1. A scalar calculation of the 
rotational momentum of the RW is given by Equation (16) [2, 
pg. 158], assuming the reference frame is body fixed to the 
RW and lies along the principal axes and the RW is fixed in 
inertial space. 
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Figure 1. Typical reaction wheel, after [4] 


h 


in 

z rw 


(16) 


Generalizing to three dimensions, while still assuming 
the same reference frame as above, yields Equation (17), 
while Equation (18) gives the most general body fixed case 
where the principal axes and reference frame axes are not 
collinear [5, p. 40]. 


- 

K\ 

V 

- 

\ 

= 

K 

i 


I x 0 0 
0 I y 0 
0 0 / 


h 

-I 


n 


I w xz 


-4 ~Iy Z 


n 


(17) 


(18) 


Applying the assumption that the RW spins only about 
the body z-axis ( a> x = co y = 0) to Equation (17) yields the same 

as Equation (1), showing that the scalar equation holds 
only when the body frame and principal axes are aligned (or 
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aligned close enough to make the products of inertia of 
negligible magnitude). 

2. Converting to the Spacecraft Body Frame 

Since Equation (17) only describes how the RW changes 
angular momentum with respect to its own reference frame, 
determining what the angular momentum change means to the 
spacecraft as a whole is the next logical step. Completing 
that step involves converting the angular momentum vector 
from the RW frame to the spacecraft body frame. Figure 2 
provides a spatial representation of the relationship 
between the two frames. 


S/C Body to RW Frame 



Figure 2. Example Spacecraft Body and RW coordinate systems 

Multiplying successive orthonormal rotations creates 
the necessary DCM. In this case, the first rotation is 
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about the 

-*™.- axis 

with a magnitude 

of 

rotation 

is about 

the 

-axis 

with 

a 

Equation ( 

;19) gives 

the 

resulting 

DCM 

b w R • 


R 


cos (-7) 

sin (-7) 

0" 

'1 

0 

0 

-sin (-7) 

cos (-7) 

0 

0 

cos(-/i) sin(-/?) 

0 

0 

1 
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-sin(-/?) cos (-/?) 
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-sin(7)sin(/?) 


sin (7) 

cos (/?) cos 

(r) 

- 

-cos(7)sin(/?) 


0 

sin(^) 



COs(/?) 



:i9; 


Using the assumptions identified above for the RW, the 
final relationship between the angular momentum in the two 
frames is given by Equation (20). 


b 1 

K 

= 

> W K 











cos(7) 

-sin(7)cos(/?) 

-sin(7)sin(/?) 


0 

0' 

X] 



= 

sin (7) 

cos (/?) cos (7) 

-cos(7)sin(/?) 
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0 
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0 

sin 

(P) 

cos(/?) 
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0 
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S^rw \ 




cos(7) 

-sin(7)cos(/?) 

-sin(7)sin(/?) 

0 




(20) 


= 

sin (7) 

cos(/?)cos(7) 

-cos(7)sin(/?) 

0 






0 

sin 

(P) 

cos(/?) 

ip. 

rw J 






-sin (7) 

sin ( P) 









= 

-cos (7) 

sin (/l) 

1 Q 

Z rw 









L cos(/?) 








The 

angular momentum in 

the body frame 

f O 

r the entire 

array 

is 

simply 

the 

summat 

ion of the angul 

ar momentum of 


each individual wheel (Equation (21)). 

n 

(2D 

i=i 

Equation (22) provides the vectorized form of Equation 


(21) . 
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~sin(;/,)sin(/?,) -sin(/ n )sin(/? (i )l 

h Krray = -cos (ft) sin (/?,) ••• - cos (/„ ) sin (/?„) : (22) 

cos (A) cos(/?„) J 1 z,F l „_ 

Given a set k of different combinations of RW wheel 
rates, an array of angular momentum points, b h array , can be 

created by a matrix form of Equation (22), leading to 
Equation (23) . 

-sin(f 1 )sin(/? 1 ) -sin(^)sin(/?,) A 0 0 

b krray = - cos (y,) si ii (/?,) ••• - cos (/„) sin (fi n ) 0 \ 0 : ••• : ( 23 ) 

cos (A) cos (A) J[o o n, h 

Use of matrix Equation (23) greatly speeds up 
execution time in the MATLAB by obviating the need for 

explicit for loops. This allows additional data points to 
be analyzed in a given amount of time, allowing for a finer 
analysis. 

For the more general case where the wheel frame is not 
aligned with the principal frame, a slight modification to 
Equation (22) is required. Although the RW still only 
spins about the z rw -axis, the coupling from the products of 
inertia, as seen in Equation (24), results in angular 
momentum being distributed along all axes in the RW frame. 
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When converting to the body frame, the entire DCM must 
now be included rather than just the third column, giving 
Equations (25) and (26). 


’h = 


b R b R ••• b R 

w rw, w rw 7 w rw„ 


X 

X 

'h 


125; 


X, = XX, 



cos(^) 

-sin(/)cos(/?) 

-sin(/)sin(/?) 

- 1 *: 

b 7 

K = 

sin(^) 

cos(/?)cos(^) 

-cos(f )sin(/7) 

X 


0 

sin(/?) 

cos(/?) 

K 


Q 


[26! 


In order to create the equivalent form seen in 
Equation (22), an augmented DCM and column array must be 
created. It will be of the form seen in Equation (27). 


’h 

array 


bn bn 

w Jv rw, w iy rw 2 


VR 

w lx rw n 


X 

X 

'h 


( 27 : 


Finally, to carry out multiple calculations 
simultaneously. Equation (27) can be expanded as shown in 
Equation (28). 


’ h = b R b R • • • b R 

array w JV rw i w iy rw 2 w iy rw n 


w 7 

K 

'h 


... X 

... "h 


"h . ••• w h 


(28; 


Creating the wheel frame angular momentum matrix can 
be done efficiently by generating it row by row (as in 
Equation (29)) vice column by column. 
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(29) 


-I 


YZ„ 


Q 


••• Q 


Note that the augmented matrices only give the actual 
momentum at a specific instant in time. Including the time 
dependency would additionally require fixing the RW frame 
to the wheel and keeping track of the angle of rotation 
over time. 


3. Angular Momentum Space 

Due to the fact that a maximum wheel speed exists for 
each RW as a mechanical limitation, there is a limit to the 
angular momentum magnitude that the array can produce in 
any given direction. Although methods exist for the 
analytical determination of the maximum momentum envelope 
[4], Equation (22) provides a basis for the rapid creation 
of a large set of momentum points within the momentum space 
based on a rigorous sweep through all combinations of RW 
speeds. A simple algorithm can then be used to find the 
maximal points for all directions around a sphere to 
determine the angular space for further analysis. 

A geometrical approach to creating this volume also 
exists [6], however the point cloud approach presented here 
allows additional flexibility when dealing with hybrid 
RW/CMG arrays. 

The momentum space itself assumes that the RW array is 
fixed to the body which is, in turn, not rotationally 
accelerating in inertial space. If the array is not fixed, 
the DCM created in Equation (19) would change over time. 
Although this can be accounted for by using the appropriate 
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time history of the angles associated with the DCM, the 
shape of the momentum space itself would change over time. 
Similarly, if the array had a rotational acceleration 
component in inertial space then this acceleration would 
have to be included in the wheel speed term, affecting the 
momentum vector direction for each RW changing the net 
momentum vector direction for the RW array as a whole, thus 
affecting the shape of the momentum space. 

Although this seems to indicate that the momentum 
envelope would, in general, not apply since the spacecraft 
is always accelerating in inertial space, the magnitude of 
spacecraft rotation is much less than the magnitude of 
wheel rotation so the effect of the spacecraft rotation on 
the momentum space of the RW array is, in fact, negligible. 

B. CONTROL MOMENT GYRO ANGULAR MOMENTUM 

1. Overview 

Although the basic principles for determining the 
angular momentum of a CMG matches those used for a reaction 
wheel, the extra degree of freedom of the CMG slightly 
complicates the equations. The rotor now rotates about the 
local x-axis, which is the CMG body-centered direction of 
the angular momentum vector, shown in Figure 3. 
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Wheel 


W z g z 



Figure 3. Typical CMG, from [4] 


The CMG gimbal rotates about the local z-axis, and 
this rotation must be accounted for when converting the 
angular momentum from the gimbal frame of the CMG to the 
spacecraft body frame. This rotation is accomplished by 
performing a z-axis rotation in the opposite direction of 
the gimbal angle, 8 , prior to doing the first rotation from 
Equation (19). The additional rotation is included in 
Equation (30). 
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cos (- 7 ) 

sin (- 7 ) 
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0 


0 

0 

1 

0 

-sin(-yff) 

cos(-/?) 

0 

0 

1 


cos(<5')cos(f)-cos(y0)sin(^')sin(f) -cos(f )sin(<5')-cos(/?)cos(^')sin(f) sin(/?)sin(f) 

cos(£)sin(f ) + cos(/?)cos (7 )sin(^) cos(y9)cos((5')cos(f)-sin(^')sin(f) -cos(/)sin(/?) 

sin(/?)sin(£) cos(^)sin(/?) cos(/?) 

(30) 


The final equation for determining the orientation of 
the angular momentum vector for a CMG in the spacecraft 
body frame can thus be determined, as shown in Equation 
(31) (or by the more general form of Equation (32)). 
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: 32 ; 


2. Angular Momentum Space 

a. Momentum Inversion Algorithm (MIA) 

Although a sweep analysis through all possible 
combinations of gimbal angles would provide a similarly 
maximal result as that performed for the reaction wheels, 
the presence of singular surfaces [7] presented an 
opportunity to develop a more unique approach that greatly 
improves the speed with which the analysis is conducted. 
Essentially, this technique is an inversion of the CMG 
array angular momentum equations typically seen in 
literature, such as Equation (33) [8, p. 466]. 
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(33) 


i =1 

-cos(/?)sin(^' 1 ) -cos(c!> 2 ) cos (/?) sin (< 53 ) cos^) 

= cos(< 5 ' 1 ) + -cos(/?)sin(<!> 2 ) + -cos(< 5 ' 3 ) + cos(/?)sin(< 5 ' 4 ) 

sin(y 0 )sin(( 5 ' 1 ) sin(/?)sin(£ 2 ) sin(/?)sin(£> 3 ) sin(/?)sin(( 5 ' 4 ) 

The basis for the MIA analysis lies in the fact 
that maximizing the amount of angular momentum each 
individual CMG places in the desired direction will yield a 
maximum result for the array as a whole. This allows the 
overall problem to be decomposed into n smaller problems, 
where the direction vector is projected onto the rotational 
plane of each CMG in the array to determine the appropriate 
gimbal angle for that CMG to maximize the momentum it 
places in the desired direction. The momentum vectors for 
each CMG are then added, resulting in a least squares 
solution to the inversion problem. Figure 4 shows the 
result of this algorithm on a four CMG array. 
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Figure 4. Schematic of the inversion algorithm 


For the i ,h CMG, consider Equation (30), rewritten 
in the standard Ax—b form as Equation (34). In this case, 

d is the unit vector along the desired angular momentum 
direction, A is the unknown gimbal angle vector, and K is 
the remainder of the matrix from Equation (30). 


cos(x) 

-sin(;r)cos(/?) 

-sin sin 

(P)~ 

cos (A) 

sin (7) 

cos(/?)cos(x) 

- cos (f) sin 

(P) 

sin (A) 

0 

sin(jff) 

cos(/?) 


0 


One problem that becomes immediately apparent is 

the fact that in Equation (34), there is only one angle 

being solved for. Using standard methods for solving the 

problem will not work since, in the absence of additional 

constraints, they will treat the two terms as decoupled and 
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solve one independently of the other. Analyzing this 
problem in a vectoral manner provides a means to move 
forward. Figure 5 represents the initial problem for a 
single CMG. 


MIA Analysis 


ro~ CMG Rotation Plane 
.Initial CMG h 

- d 



Figure 5. Gimbal rotation plane and desired direction 

vector in the CMG frame 


Geometrically, the left hand side of Equation 
(34) describes a plane in the spacecraft body frame. 
Figure 5 provides a visual example of this, with the circle 


describing the 

possible 

orientations 

of 

the CMG 

angular 

momentum vector 

and d 

describing 

the 

desired 

angular 


momentum direction. If d lies in the rotation plane, the 
CMG can orient itself such that its angular momentum vector 

lies parallel to the direction of d . If d does not lie in 
the rotational plane, then the ideal direction for the 

angular momentum vector would be along the projection of d 


25 














onto the plane. To do the projection, the normal ( n + ) and 
anti-normal ( n ~) vectors are calculated based on the 
momentum directions at two orthogonal gimbal angles, such 

as 0° and 90° ( A( 0 °) and A( 90 °) respectively), shown in 

Equation (35) and diagrammatically by Figure 6. 


n + = ^A( 0 °)x^A( 90 °) 
n- =KA(90°)xKA(0°) 


(35) 


< t 


MIA Analysis 





CMG Rotation Plane 


- CMG h 

- CMG h 


fro 
5=9 o' 


Figure 6. Finding the Rotation Plane Normal Vectors 

The appropriate normal or anti-normal vector is 
chosen such that the relationship given in equation (36) is 
satisfied. 

h*d >0 ( 36 ) 

The projection can now be determined by 
subtracting the perpendicular component in Equation (37) 
and Figure 7. 
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(37) 



< \ 


MIA Analysis 
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d 

" d projection 
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i 


Figure 7. Projecting d onto the Rotation Plane 


The least squares gimbal angle can thus be 
calculated using Equation (38), as described in Figure 8. 


8f 


cos -1 (^TA (0 o ) m d proj ) if(KA{0°)xd proj )-n + >0 
360° -cos -1 ^KA(0°yd proj j otherwise 


(38) 
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MIA Analysis 



ro~ CMG Rotation Plane 
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Figure 8. Determining the Gimbal Angle of the Projection 

The resulting CMG momentum vector described by 
placing the CMG at the calculated gimbal angle provides the 
least squares solution for that particular CMG. Orienting 
the CMG at any other gimbal angle would increase the norm 
of the difference between the normalized CMG momentum 
vector and the normalized direction vector. Performing the 
same algorithm on each of the other CMGs in the array 
provides the individual least squares solution for those 
CMGs as well. Summing the momentum vectors provided by 
those solutions gives a net momentum vector for the array. 

As seen in Figure 4 there may be some angular 
difference (or angle residual) between the desired 
direction and the least squares net momentum vector. It is 
possible to construct equations that calculate this angle 
residual as a function of solved gimbal angles. This 
equation can be used with a numeric solver to reduce the 
residual to near zero and thus determine a quasi-least 
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squares gimbal angle solution for any given direction in 
which a solution exists (a solution will not exist, for 
example, in a direction that is perpendicular to the 
rotational planes of all CMGs). In effect, the preceding 
algorithm simply provides an initial guess so that a 
numeric solver can find a set of gimbal angles that give a 
near maximum net momentum vector for a CMG array. This, 
therefore, allows for the construction of a quasi-maximal 
surface by sweeping around a sphere rather than sweeping 
through all possible gimbal angle combinations. 


b. MIA Angle Correction 

Assuming unit momentum, the momentum vector for 
each CMG is calculated via Equation (39), using the gimbal 
angle given by Equation (38). 

cos(y) -sin(^)cos(/?) -sin(^)sin(/?) cos(^) 
sin(^) cos(y9)cos(^) -cos(/)sin(/?) sin(j) = b h cmg (39) 

0 sin(/?) cos (/?) 0 


Since all vectors are known at this point, the 
angle residual can be determined using the dot product 
definition in Equation (40) . 



(40) 


If 0 res falls outside of a given threshold, the dot 
product itself can be used to solve for the gimbal angles 
that reduce the residual to near zero based on the identity 
given in Equation (41) . 


u || v 

A A < 

u»v = 1 


(41) 
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Therefore, when the angle residual is zero. 
Equation (42) will be true. 

d» b h -1 = 0 (42) 

Rapidly finding the solution to Equation (42) is 
well within the capability of modern numerical solvers, 
assuming the solution exists. 

If the products of inertia are non-zero for the 
CMG wheel, the angular momentum vector will be offset from 
the spin axis of the rotor. This will induce an error not 


taken into account by 

the 

equations 

as 

developed. In 

general, 

the products 

of 

inertia 

will 

be orders of 

magnitude 

smaller than 

the 

moments 

of 

inertia, greatly 


reducing the size of the error. To correct for this error, 
the projection of the angular momentum onto the rotational 
plane should be used in place of K A . 

c. Momentum Space Analysis 

A maximum momentum surface for the CMG array 
based on the MIA can now be quickly constructed. Similar 
to finding the maximum surface for RW arrays given a point 
cloud, the first step involves creating a unit sphere with 
the point density determining the fidelity of the CMG 
surface construction. Iteration through the spherical 
envelope provides the direction vectors for the algorithm. 

The output of the algorithm provides both a 
maximal surface for all analyzed directions, as well as an 
angle residual which gives a measurement for how well the 
CMG array can point its momentum vector in a given 
direction. 
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C. SUMMARY 

This chapter developed the equations and algorithms 


necessary 

to create maximum momentum volumes for 

both 

RW 

and CMG 

arrays. 

These 

momentum volumes can be 

used 

to 

estimate 

the 

maximum 

theoretical performance 

of 

an 

arbitrary 

array 

and for 

more advanced techniques 

such 

as 

the analysis of 

hybrid 

arrays and performance comparison 

between 

gimbal 

angle-limited singularity free CMG array 

momentum 

volumes 
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IV. TORQUE SPACE 


A. OVERVIEW 

The general equation of motion for a spacecraft is 
given by Equations (43) and (44) [8, p. 463], where J is 

the total moment of inertia of the system, CO is the 
spacecraft body rate vector and b h MED is the MED momentum 
vector expressed in the spacecraft body frame. 

b H s +co b,N x b H s = b T ext (43) 

b H s = Jco b,N + b h MED (44) 


Substituting Equation (44) into Equation (43) and 
assuming a torque-free environment yields Equation (45) [8, 

p. 463]: 

( Jd> b,N + b h MED j + co blN x(jco hlN + b h MED ) = 0 . (45) 

Equation (45) can then be rewritten to separate the 
spacecraft terms and the MED terms, shown in Equations (46) 
-(48) . 




b/N I „ B/N s, ls-P' N 
+ 00 xj CO 


+ 

Js/c L 


b U + m blN x b h 

,l MED X ,l MED 


= 0 


J MED 


' 46 ; 


b h = b h 

MED n 


+ b h 


(47) 


b h = b h + b h (48) 

Recalling Equations (16)-(18) and applying them to 
Equation (47) yields Equation (49). Equations (46)-(49) 
hold regardless of the reference frame they are expressed 
in, so long as all terms are expressed in a consistent 
frame. 
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^MED Irvfc^rw + ^cmg^r ( 49 ) 

Equation (49) allows Equation (46) to be further 
qeneralized into Equation (50). 


~Jeol 


■ b,N +co b,N xJco b,N ~ 


S/C 


+ 


I co n "' N + /, (b rwlN +1. co cmg,N +1 cb cmg,N + (0 h,N X 


cmg 


cmg 


(/ co IWIN +I 

0 J cnglN ) 

\ rw 

cmg J 


:5o: 


= o 


MED 


Equation 

(50) can be 

expressed 

by Equation 

(51) 

by 

assuming that 

I RW =0 in the 

spacecraft 

body frame. 

This 

is 


a reasonable assumption as reaction wheels will be nearly 
axisymmetric about their spin axis, which is fixed in the 
body frame. Note that there is also the implicit 

assumption that the overall inertia of the spacecraft, J , 
is much qreater than I CMG r so it will remain rouqhly 
constant. Since the CMG terms are much more complex, they 
will be folded back into quantities that will be developed 
separately [8, p. 464]. 


" J(b b,N +co blN 


, h cmg +a b,N x b h cmg =- b r cmg 


1 

r t +t i 

Js/C 

rw cmg J 


= ~ b T„ 


(51) 


B. REACTION WHEELS 
1. Derivation 

Althouqh Equation (51) provides a method for 
determining the inertial torque for the MED arrays, 
analyzing the torque capability in the body frame (as was 
done for momentum in Equation (23)) removes the 
complication of having to include inertial motion. 
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Equation (52) provides the relationship between momentum 


and torque within the local 

wheel frame. 





II 

r 'l 


(52) 

Recalling 

Equation 

(16), Equation 

(52) 

becomes 

Equation (53): 

tv _ W 1 T 

r = L 7 ™ 

„ rw/N , t •rw/N~ 1 

co + I nv co J . 


(53) 


The RW is axisymmetric about its spin axis, so the 
wheel inertia will remain constant. Converting to the body 
frame is done in a similar matter as done for momentum, as 
given by Equation (54) . 

b r= b w R ] + HR w [/ w © w/w ] (54) 

Since the wheel frame is fixed with respect to the 
body frame, the DCM remains constant, yielding the final 
torque equation. Equation (55), for a single RW. 

b T=ZR w [l rw j> rw,N ] (55) 

Assuming the body is not accelerating in inertial 
space. Equation (55) can be reduced to Equation (56), 
which, as expected, is the time derivative of Equation (20) 


b 


r = 




0 

0 

Q 

rw 


(56) 


To determine the torque generated by the entire array. 
Equations (21) and (23) can be differentiated and rewritten 
as Equations (57) and (58) . 


35 




-sin(f 1 )sin(/? 1 ) -sin(/ n )sin(/? H ) I Zl 0 0 Q„ ^ 

Murray = -cos (>, )sin (/?,) ••• - cos (/„) sin (/?,) 0 0 : - : (58) 

cos(^) cos(/? n ) J[o 0 I z J Q ;ii 

A maximal torque space can then be developed via the 
same method used to develop the maximal momentum space. 

Since the geometry of the RW array does not change 
with respect to the body, the shape resulting from Equation 
(57) will be similar to the shape developed by Equation 

(23) . 

C. CMGS 

Although the torque-momentum relationship given in 
Equation (52) holds for CMGs as well, the resulting 

derivative is more complex due to the time variant gimbal 
angle, 8 . This additional degree of freedom causes the 
available torque to become state dependent. In other 
words, each individual CMG configuration has its own torque 
space. Among other problems, this makes visualization 

difficult as the torque currently available to the CMG 

array is not a constant function of the gimbal angle and 
will change with time and for each gimbal trajectory. 

Traditionally, this issue causes most analyses to be 
conducted on the momentum space of CMG arrays, and 
significant work has been done in that area [ 9; 10; 11; 12; 8, 
pp. 466-470] . In particular, much of the work has focused 

on singularity analysis and either developing control laws 
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to avoid singularities or to create CMG arrays that 
minimize singularities [12]. An alternative analysis 
involves looking at the state space of the array. 

Developing a state space analysis for CMG arrays is 
numerically intensive, especially if a high fidelity is 
desired. Furthermore, there is an issue of dimensionality. 
By definition, a state dependent plot of available torque 
would describe the momentum output of the array for each 
gimbal angle and gimbal rate. When considering a four 
wheel array, this implies there are eight input dimensions 
that must be mapped onto three output dimensions. 

Despite the problems associated with working in state 
space, there are several benefits that can be realized. 
For example, by identifying regions of the space that are 
composed exclusively of non-singular states it is possible 
to constrain the allowable gimbal angles for the array so 
that the basic pseudoinverse control law is sufficient to 
maneuver the satellite. Without constraining the gimbal 
angles, this approach would normally result in the loss of 
control of the spacecraft after a short period of time due 
to singularities. Another option would be to select a 
larger region by allowing only hyperbolic singularities to 
be included, as these can be avoided by standard techniques 
such as adding null motion to the pseudoinverse control law 
[ 8, p. 4 68]. 

1. State Space Analysis 

Equation (59) represents one of the key assumptions 

for this analysis. This assumption allows the input 

dimensionality to be reduced by half, greatly shrinking the 

scope of the problem. There is still the issue of 
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visualizing data that lives in R" +3 , however, where n is the 
number of CMGs in the array. By constraining the output of 
the analysis to be a scalar rather than a vector, the three 
output dimensions can be further reduced to one. This will 
result in the loss of some information; however the 
information may be unnecessary, depending on the particular 
analysis performed. 

5^5 (59) 

2. Contiguousness 

At this point, it is necessary to introduce the 
concept of contiguousness. In this thesis, a state space 
analysis representation is considered contiguous if all 
continuous trajectories that can be traced in the state 
space are physically realizable. When analyzing a single 
CMG, it would require some effort to present a state space 
visualization that is not contiguous. When reducing 
dimensionality during a CMG array analysis, however, it is 
very important to keep this property in mind. When 
considering the problem at hand, developing gimbal angle 
limits to achieve a singularity free operational momentum 
space for a CMG array, a non-contiguous analysis may allow 
for a state within this space that can only physically be 
reached by violating gimbal angle limits. 

For example, consider a four CMG array. Using the 
previously discussed assumptions, the state space of this 
array can be reduced from R" to R~" . To convert this down 
to R 3 , the gimbal angles for two of the four CMGs can be 
represented by an index number along the X and Y axes of a 
topographic plot. This state index number refers to a 
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specific state combination of the two CMGs chosen for that 
axis. Assuming that adjacent state index numbers represent 
gimbal angle combinations that are physically next to each 
other, the resulting plot will give a semi-contiguous 
representation of the system, shown in Figure 9 and Figure 
10. The remaining two CMGs in the array would be 
represented in state index two and have a similar plot to 
Figure 10 . 


State Space Analysis Sample Plot 



State Index 1 


Figure 9. Four CMG semi-adjacent state space 
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State Index 1 States 



Figure 10. Reduced dimensional state space analysis state 

index array with two CMGs 

This semi-cont iguousness arises from the loss of one 
degree of freedom in each axis. This implies that although 
any trajectories drawn in this representation will result 
in valid physical trajectories for the CMG array, it does 
not represent the full set of possible trajectories. To 
achieve this, each CMG has to be able to change 
independently, requiring a plot of greater than three 
dimensions, which is very difficult to visualize. 

Aside from indicating the problem with dimensional 
reduction and maintaining full contiguousness. Figure 9 and 
Figure 10 provide a guide for how to determine gimbal angle 
limits. Assume that the plot in Figure 9 shows locations 
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that result in singular (black) and non-singular (white) 
states in the CMG array state space. Gimbal angle limits 
for singularity free space can now be pulled from that plot 
by finding state index limits that bound only white areas. 
Specifically, if the gimbals were limited to those 
represented by the gimbal angles in states one and two for 
state index one and five through eight for state index two, 
the resulting momentum space for the array would be 
singularity free. By Figure 10, the first CMG in the 
array, as defined by the user, must operate between zero 
and one radians while the second CMG must remain at one 
radian. An additional plot for state index two would be 
required to determine the gimbal angle limits for the rest 
of the CMGs. 

D. SUMMARY 

This chapter developed the equations necessary to 
determine the torque envelope for a RW array and showed 
that the geometry of the torque envelope will have the same 
shape as the momentum envelope. 

For CMG arrays, the groundwork was laid to use a state 
space based analysis to examine the possibility of 
operating in singularity free momentum space. The concept 
of contiguousness was established to ensure that the 
results of the analysis would allow for the CMGs to reach 
every state in the defined non-singular region without 
having to violate the established boundaries. Subsequent 
chapters will look at how the limited momentum envelope 
compares to the theoretical maximum developed in the 
Chapter III, and test the feasibility of operating within 
this region. 
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V. 


ANALYSIS OF REACTION FORCES 


One important area that has received little attention 
in the literature is examining the reaction forces at the 
attachment points of the RWs and CMGs throughout a 
spacecraft maneuver. Current pedagogy on this topic 
focuses on the transfer of momentum between the spacecraft 
and the MED. The vehicle of this transfer, reaction forces 
and moments on the attachment points between the MED and 
the spacecraft, is rarely discussed. There are several 
possible causes for this. One reason may be that 
experimental data determined that standard attachment 
design was robust enough. Or that other loading such as 
vibrational loading greatly exceeded the expected magnitude 
of maneuver loading. It is even possible that the analysis 
sits in a particularly bad convergence between two fields 
such that neither dynamicists nor structural engineers have 
a desire to study the problem. Regardless of the reason, 
it is important to understand what is happening at the 
attachment points in order to appropriately size the 
motors, bearings and gimbals so that total mass can be 
minimized while maintaining an appropriate safety factor to 
protect against failure. 

For the analysis presented in this chapter, all motion 
between adjacent surfaces is assumed to be frictionless and 
all forces on bearing surfaces are assumed to be acting as 
a single resultant force at the centroid. The effect of 
gravity is neglected. All final reactions are defined as 
acting on the component being analyzed. 
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A. 


INERTIAL MOTION 


Solving for the reaction forces for the MEDs begins 
with determining the motion of each component part within 
inertial space. To calculate the motion specific to these 
components, data such as satellite orientation, gimbal 
angles, etc... must be known. The rates and accelerations 
associated with these motions can be estimated using 
numerical differentiation techniques, but the fidelity of 
the analysis may suffer. The most accurate results require 
complete knowledge of the system during the time period 
analyzed. 

The inertial motion for both the RW and the gimbal 
frame of the CMGs is the same. To begin, the inertial 
position vector for the RW center of mass must be expressed 
as in Equation (60) . 


N 


v r) b 

Ps/c r ™- 


(60) 


As noted in the introduction, there are two equivalent 
methods to proceed from here. The first method, shown in 
Equations (61) and (62), uses the transport theorem and 
differentiates twice, then uses the appropriate DCMs to 
convert all terms to the correct frames. 


W v„=^c+(/ W xV s/c ) + t + rxV„, 
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'a = N P slr + h r +2(co blN x h r ) 

rw S/C rw \ rw J 


OJ h,N x V , + a hlN x V , + co blN x co blN x V 


To maintain all vectors in a consistent basis, all 
body terms must be expressed in the inertial frame. 
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Alternatively, Equation (60) can be expressed so that 
all quantities are expressed in the inertial frame and the 
differentiation operations performed (Equations (63) and 
(64) ) . 
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1 SIC ^ b 1X 'rw 
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Then, using the identities in Equation (65), Equation 
(64) can be written in the same form as Equation (62) as 
shown in Equation (66) . 
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Equation (66) gives the same expression as Equation 
(62) . The primary difference between the two approaches 
lies in the fact that the product-rule approach, the second 
approach, requires that the vector be resolved in the 
inertial frame prior to differentiation while the transport 
theorem approach does not. This can cause the product-rule 
approach to create a large amount of terms as additional 
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frames are added to the problem, making it fairly unwieldy 
to solve by hand. The primary benefits of the product-rule 
approach deal with the simplicity in creating an algorithm 
to do symbolic differentiation on a computer as well as 
explicitly maintaining all terms in a consistent reference 
frame. In this case, there are sufficiently few reference 
frames that both approaches are reasonable to work with. 


If the RW center of mass is assumed to be fixed within 
the structure, a valid assumption for a rigid spacecraft, 
then some of the terms can be eliminated to give: 


'^= n p s ,c+:r 


b/N 1 . r „b/N 1 r „b/N 

a +1 eg eg 
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As mentioned earlier, the CMG gimbal frame has the 
same equations of motion as the RW frame, rewritten as 
Equation ( 68) : 


N 
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r z 
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= P + Rr 

r S/C ^ b ,X 'g 


= N Ps,c + N b R 


r nb/Nl . V 

a + (o (o 



( 68 ) 


For a CMG, the presence of an additional rotor frame 
adds another term to the initial position vector and 
requires an additional DCM to go from the gimbal frame to 
body frame. Using a similar analysis as above and applying 
the rigid body assumptions yields the rotor equations of 
motion (Equation (69)). 
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N a r = N a g + N b R{[q hlN ] + [co blN ][cg blN ]) h g R g r r 

+2 N b R[d’ ,N ] b g R[co 8 ' b ] g r r +2 N b R[<p b,N ] b g R 8 r r 

+2^;/?[ffl* /6 ]*r r + ^^([g i/i ] + [®* /6 ][®* /6 ])*r r (69) 

+ N b R b g R g f r 

N a r = N a g + N b R([q b/N ] + [^ N Jcu b/N ]) b g R g r r 

+ 2 N b R[a b,N ]*R[g g,b ] g r r + N b R>R([q g,b ] + [a* ,b ][g s,b ]) s r r 


When using the transport theorem to develop Equation 
(69), an additional step is required to determine the 
angular rates, since it is very unlikely that the inertial 
angular rate and acceleration of the gimbal frame is 
measured. This quantity can, however, be determined from 
the spacecraft body rates and the gimbal rate (shown in 
Equation (70)), both of which are typically measured 
quantities . 
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B. REACTION WHEELS 

1. Case One - Cantilever 

For RWs, two cases have been considered, as 
represented by Figure 11 and Figure 12. For case one 
(which represents a cantilever arrangement), there is 
exactly one reaction force and one reaction moment, leading 
to Equations (71) and (72) . 
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Figure 11. Schematic of RW cantilever arrangement 
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(72) 


The RW motor, assumed to be mounted with bearing one, 
is not explicitly included on the left hand side of 
Equation (72). It is, however, implicitly included in the 
^inertial term as the effect of that torque on the system: d> RWz . 

If it is desired to use actual motor torque output instead 
of measured or calculated wheel acceleration, the 
additional torque should be included. The wheel 
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acceleration due to the motor torque, however, should then 
not be included in the equations of motion. Doing so will 
cause the terms to cancel out and remove the motor torque 
from the overall reactions. 

The right hand side of the force balance equation is 
straightforward as the acceleration term was developed in 
Equation (67) and all quantities in that equation are 
known. The right hand side of the moment balance equation 
requires some additional analysis to be put into a useful 
form. 

Since the rotation of the RW frame in inertial space 
is not directly measured, it must be calculated from known 
quantities as shown by Equation (73) . 

rw/ b/ rw/ 

(j/N = ™ Rco /N + 0 j/h (73) 


As noted in Section A, to utilize the product-rule 
approach of differentiation it is important to initially 
express the vector in the inertial frame. After 
differentiation, the frame may be changed as desired. 


The angular momentum of the RW can now be calculated 
and differentiated (Equations (74) and (75)) . The moment 
of inertia for the RW is assumed to be constant about the 
spin axis . 
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The third term of Equation (75) includes the matrix 
^<y Ww J. Physically, this matrix is the rotation of the body 


with respect to the RW frame, given in body frame 
coordinates. Since it is a skew symmetric matrix, there 
are two equivalent approaches: change the basis of the 
rotation vector prior to converting it to a skew symmetric 
matrix, or use a similarity transformation to transform it. 
Since the data will be given in the form of a rotation rate 
in the RW frame, the first method should be used. Thus, 
Equation (76) now contains terms which designate given or 
easily calculated data. 
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Finally, converting into the RW frame yields Equation 

(77) . 

™h = r "R b R N h 
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The motor torque is equal to the third component of 

Equation (77) [1, pp. 179-185]. Under ideal conditions, 

the motor would resist any change in wheel speed with 

respect to the body due to the motion of the spacecraft 

around the RW spin axis unless a torque input is desired. 
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The inclusion of an additional motor torque term in the 
left hand side of the torque balance equation would 
(assuming ideal data) negate the spin-axis torque and 
corrupt the output. 

2. Case Two - Fixed-Fixed 

The solution for case two (given in Figure 12) 
requires some additional assumptions due to the 
overdetermined nature of the fixed-fixed system. These 
assumptions, along with solutions based on continuity and 
compatibility conditions allow for determining the bearing 
reactions. 


Axial 



Figure 12. Schematic of RW fixed-fixed arrangement 

For forces along the axial direction, it is assumed 
that the non-motor bearing, bearing two, is allowed to 
translate and thus cannot support a load. This places all 
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axial loading on bearing one. Similarly, for torques in 
the axial direction, bearing two is assumed frictionless so 
all axial reaction moments are resisted by the motor at 
bearing one. 

For the transverse cases, reaction forces and moments 
for each axis can be calculated via superposition using the 
solved cases in a reference book such as Table 8.1 of 
Roark's Formulas for Stress and Strain [13, pp. 208-223]. 


To 
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142] . 


Additionally, all forces and moments must be resolved into 
the RW frame prior to the analysis. 
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The final reactions are given in Equations (80)-(85). 
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C. CMGS 

In order to simplify the analysis of a real CMG, such 
as the one shown in Figure 13, a model was constructed with 
a rotor shaft mounted directly onto a gimbal shaft (Figure 
14). The only other assumption made in this simplification 
involved the end conditions of both shafts. All other 
characteristics such as length, center of mass location, 
mass and moment of inertia were left as variable. Neither 
the actual rotor nor gimbal frames need be in their 
respective principal frames. Instead, those body frames 
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are defined with respect to the geometry of the system. If 
the body frame and principal frame are not aligned for the 
rotor, the system is assumed to be static. 


Spin Axis 
Reference 



Gimbal 

Axis 


Figure 13. Typical CMG Construction, from [14] 
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Axial 



BRG1 


BRG2 


Gimbal 




Rotor 


BRGl 


BRG2 




BRGl 

CMG 


Figure 14. Simplified CMG model 


The CMG analysis is broken up into two parts: analysis 
of the rotor and analysis of the gimbal. The analysis for 
the rotor matches the previous RW analysis very closely. 
The gimbal, however, adds additional complexity due to the 
mounting of the rotor on the gimbal shaft in this model. 
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Both the rotor and gimbal are independently analyzed 
as a cantilever and as a fixed-fixed system, giving four 
possible configurations for the CMG as a whole. 

1. Rotor 

a. Cantilever 

Equations (71) and (72) both apply to the 
cantilever CMG rotor (shown in Figure 15), with the 
inertial acceleration already determined by Equation (69). 
Equation (86) gives those rewritten equations. 

F r = m r a T 

M R A + rX { F ,ner,ial ) = 1 ^ X K ^ 

t? ( 86 ) 

F Inertial = m r a r 

M Inenial = ( 7 ,. ^ +0)‘ /n X I r a/ N ) 


Figure 15. 




◄-► 

t 


Schematic of CMG rotor cantilever arrangement 


All that remains, then, is to derive the 
derivative of the inertial angular momentum vector to 
create the inertial moment. The angular momentum vector 
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itself has already been discussed and is given in Equation 
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Expanding the angular rate term and writing the 
vector in the inertial frame yields Equation (88). 


A h r = b R b R g r R /, :;R g b Rof' N + N b R b R g r RI r 'R co g ' b 

( 88 ) 

+ N b R b R g r RI r co r,g 

Further simplifications occur on the realization 
that g r R-' g R-I 3x3 and result in Equation (89). 
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Differentiating as above and converting back to 
the rotor frame gives Equation (90). 
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b. Fixed-Fixed 

The same force and moment balance equations used 
by the fixed-fixed RW (Equations (80)-(85)) also apply to 
the rotor in Figure 16. , however the transverse directions 

are now the y-z plane and the rotation axis is now the x- 
axis, so only a slight modification of the subscripts is 
required, which are shown in Equations (91)-(95). 
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Figure 16. Schematic of CMG rotor fixed-fixed arrangement 
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2 . 


Gimbal 


a. Cantilever 

For the gimbal. Figure 17, an additional force 
and moment must be considered: the reaction of the rotor on 
the gimbal. As these quantities are known from the rotor 
calculations, they simply add additional terms when solving 
for the cantilever reactions, given in Equation (95). 




Figure 17. Schematic of CMG gimbal cantilever arrangement 
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jb. Fixed-Fixed 

To compute the gimbal reactions for the fixed- 
fixed case (Figure 18), the rotor force and moment must be 
shifted to the gimbal''s center of mass. Equation (96), in 
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order to satisfy the condition that only a single force and 
moment are applied to the beam. 



Figure 18. 
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Using the idea of the combined force and moment, 
the final solutions are nearly identical to those presented 
in the fixed-fixed RW case and are given in Equations (97)— 
( 101 ) . 


= F* 

Inertial z ’ 

M r 

*BRGlz 

= M* 

Inertial z 

= 0; 

M r 

*-BRG2z 

= 0 


(97) 


Inertial brq\ x ^ 


-( l-a)\l + 2a ) 


r Inertial brg 2 


°^a 2 (3l-2a) 


( 98 ) 
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D. SUMMARY 

This chapter developed the equations of motion and 
applied forces and moments necessary to analyze the 
reaction forces of MEDs of various endpoint configurations 
over the course of a spacecraft maneuver. The derivation 
of these equations focused on requiring only the telemetry, 
geometry and operating characteristics that would be 
readily available from either an actual satellite or 
simulation. 

The validation of the developed equations will occur 
in Chapter VI, Section C. This validation will show that 
the developed equations behave properly under various 
assumed condition. Following validation, simulated 
spacecraft maneuvers will be analyzed, with the focus on 
showing the impact that different control laws can have on 
the reaction forces during similar maneuvers. 
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VI. CASE STUDIES 


In this chapter, a number of case studies are used to 
show the results of analyzing standard sets of both RW and 
CMG arrays using the algorithms developed in the previous 
chapters. Additionally, the performance of non-standard 
analyses such as sensitivity analyses on RW balancing with 
respect to the array momentum envelope are included to show 
the flexibility of the algorithms. 

A. MOMENTUM ENVELOPE 

1. Reaction Wheels 

For the following RW analysis, the wheel 
characteristics are based off of the Bradford Engineering 
W18 RW [15] . The important characteristics of the RW are 
as follows: 

co RW = + / - 4000 RPM =419 rad / sec 
d) RWi = + /- 55 rad /sec 2 
m -5 kg 

".0022 0 0 
I RW = 0 .0022 0 kg ■ m 1 

0 0 .0045 

r RWmia = ± 55(.0045) = .2475 N-m 

a. Nominal Momentum Space 

When an RW array is composed of a single RW, the 
available momentum space forms a line that lies along the 
spin axis, as shown by Figure 19. 


63 



One RW Momentum Space 



Figure 19. Momentum space for a single RW 

For three RWs in an orthogonal configuration, the 
momentum space forms a box where each of the eight vertices 
is obtained by operating, in one of eight combinations, 
each RW at its fastest speed in either the forward or 
reverse direction. Figure 20 shows an example of this 
shape. 
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Three RW Momentum Space 



Figure 20. Momentum space for an array of three RWs in an 

orthogonal configuration 

Adding a fourth RW and orienting all four RWs in 
the standard 54.73° pyramidal orientation [6], with y and 
P angles given by Table 1, creates a more complex shape: a 
dodecahedron. 


RW 

Y ( de 8) 

P ( de g) 

1 

0 

54.73 

2 

90 

54.73 

3 

180 

54.73 

4 

270 

54.73 


Table 1. Parameters for a four RW 54.73° pyramidal 

configuration 
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As before, the vertices represent different 
combinations of wheels operating at their maximum speed in 
either the forward or reverse direction. The dodecahedron 
is shown in Figure 21. 


Four RW Momentum Space 



Figure 21. Momentum space for an array of four RWs in a 

54.73° pyramidal configuration 

Although there is more than one configuration of 
RW arrays, the characteristic of a given array will 
ultimately dictate its suitability for a given application. 
For the three RW box array, each dimension is driven by a 
single RW, so the size of each side is equal to twice the 
momentum capability of each wheel. For the four RW 
pyramidal array, the skewed RWs allow additional RWs to 
affect the array momentum, even if the desired direction is 
along the rotation axis of one of the RWs. Ultimately, 
this means the largest sphere that can fit inside the 

momentum volume for the pyramidal array is larger than that 
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for the box array. The theoretical maximum performing 
array according to the maximum sphere metric would be one 
that has a perfectly spherical momentum space. This would, 
however, be impractical to build as it would require a 
large number of devices. 


b. Skewed Momentum Space 

One benefit to calculating the momentum space 
using Equation (28) instead of Equation (22) is that wheel 
imbalances can now be taken into account. 


Assuming that the same RW as analyzed above now 
contains a large imbalance such that the body z-axis of the 
wheel is offset from the principal z-axis by 10°, the 
single RW momentum space line becomes skewed as shown by 
Figure 22 when frozen at a particular instant in time. As 
the wheel rotates, the momentum vector will create a cone, 
an effect commonly called nutation [5, pg.53]. This makes 
the momentum vector time dependent and complicates the 
analysis. 

One possible skewed inertia matrix is 


2.235rl0 3 -3.468.xKT 5 

-3.468.xKT 5 2.235.xKT 3 

2.78 IjcIO 4 -2.78 LxKT 4 


2.78 l.xKT 4 


-2.78 LxKT 4 
4.431x10 3 


kg - nr . 
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One RW Momentum Space (Nominal) One RW Momentum Space (Skewed) 




Figure 22. Balanced (left) and unbalanced (right) RW 

comparison 


When applied to a three or four wheel array, 
these unbalanced wheels skew the entire momentum envelope. 
This is clearly shown in Figure 23. These shapes only 
apply when the noted requirement for the applicability of 
Equation (28) is met: that the shape is valid only for a 
particular instant in time. If examined over time, the 
shape would "breath" as each vertex moved with the momentum 
cone for each RW causing the faces to change size, shape 
and orientation. 

Although the 10° offset makes the differences in 
the envelopes clear, the offset represents an imbalance 
that would result in a shortened component lifetime. More 
realistically, the offset would likely be much less than 
1 ° . 
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Three RW Momentum Space (Nominal) Three RW Momentum Space (Skewed) 



h y (N-m-sec) 44 h x (N-m-sec) h y (N-m-sec) 44 h x (N-m-sec) 

Four RW Momentum Space (Nominal) Four RW Momentum Space (Skewed) 



Figure 23. Balanced (left) and unbalanced (right) multi-RW 

array comparison 

One question that arises from this study is what 
effect this skewed shape has on the total enclosed volume. 
A smaller enclosed volume represents a performance loss, so 
determining the degree of this loss in the range of 
expected imbalance represents a useful exercise. By 
iterating the offset angle and calculating the resulting 
momentum envelope volume, a relationship can be 

established, as shown in Figure 24 and Figure 25. 
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Three RW Array Imbalance Analysis 



Figure 24. Three RW array imbalance sensitivity analysis 


Figure 25. 


Four RW Array Imbalance Analysis 



analysis 


The results show less than a .04% performance 
drop at 1° skew. The small reduction in performance for 
what is still a large imbalance indicates that assuming the 
RW rotates about a principal axis will have a negligible 
impact on preliminary ACS design. 
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2. 


CMG Momentum Space 


Consider the RW analyzed in the previous section now 
mounted on a gimbal. Further, assume that the wheel only 
rotates at its maximum speed in the forward direction and 
that the axis of rotation is the rotor x-axis, to be 
consistent with the system developed earlier. The 

available momentum space for a single CMG is a circle, as 
shown in Figure 27. 

Figure 28 shows the angle residuals between the 
desired momentum vector and the momentum vector that the 
CMG can produce via the MIA process described in Chapter 
III, Section III.B.2. The desired direction vector is 
determined from the azimuth and elevation angles around the 
unit sphere. Figure 26 provides the convention to 

determine the direction of d using the (j) (elevation) and 9 
(azimuth) angles. It is clear that this single CMG array 
can only match the desired vector when it lies in the plane 
of rotation. In other words, as the elevation component of 

d deviates from 0°, the residual from the MIA produced 
momentum vector will be equal to that deviation. 
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4> and e Convention 



Figure 26. Convention for using (j) and 0 to describe d 


One CMG Momentum Space 



Figure 27. Momentum space for an array of one CMG 
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Angle Residuals 



S(deg) 0 0 * (deg) 


Figure 28. Angle residuals between the desired momentum 
vector and the actual momentum vector 

Placing three CMGs in an orthogonal configuration 
gives a much more interesting shape, seen in Figure 29. 
Note that this figure represents the maximum momentum 
theoretically achievable in a given direction, however due 
to the presence of internal singularities it may not be 
physically possible to reach the surface in a given 
maneuver. 
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Three CMG Momentum Space 



Although a rigorous singularity analysis was not 
performed, it is likely that the conical shapes match the 
shape of a singularity surface [4] due to the relatively 
large angle residual shown in the vicinity of 
Figure 30. 


Figure 29. Three CMG array maximum momentum space. 


h y (N-m-sec) 

Three CMG Momentum Space 
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Three CMG Momentum Space 


4 - 3 - 2-101234 

h x (N-m-sec) 


-2 0 2 

h x (N-m-sec) 


each by 
















Angle Residuals (deg) 



<(> (deg) 



Figure 30. Three CMG array angle residual plot 


For the 54.73° CMG pyramid, the standard shape appears 
in Figure 31. The large vertical stripes present at 0=90° 
and (f> = 210° in Figure 32 are likely due to a singularity 


present when 



In these directions, changing 6 has 


no effect on the direction of d . This means that as the 
MIA algorithm sweeps through the azimuth values, the same 
residual will be seen. Further analysis is required to 
determine why the vertical stripe is not preset at 0 = 90° 
for Figure 30. 
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Four CMG Momentum Space 



Four CMG Momentum Space 


Four CMG Momentum Space 



Figure 31. 


Four CMG array momentum space 
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Figure 32. Four CMG array angle residuals 


3. Hybrid RW/CMG Array 

By maintaining all momentum data in point clouds, it 
becomes possible to analyze the maximum momentum surface of 
a composite array with both RWs and CMGs. To the best of 
the author's knowledge, such an analysis has not been 
reported in the literature. Figure 33 shows such an array, 
with three CMGs and three RWs oriented in two sets of 
orthogonal box arrays, similar to the three device cases 
analyzed previously.. There is some distortion of the CMG 
shape (beyond the increased size from adding the volumes). 
Internally, the CMG singularity surfaces get pushed to the 
border of the original RW surface, as that provides a 
singularity free volume. 


77 











This singularity free space implies that it may be 
possible to include a separate RW array on a spacecraft for 
the purpose of elliptic singularity avoidance. The RWs 
would allow a torque-generating CMG motion to bypass the 
singularity threshold, potentially allowing for an expanded 
utilization of the CMG performance envelope, or the ability 
to use a simplified steering law. 

3 CMG/3 RW Combined Momentum Space 


< 1 > 

</> 

E 

z 



hy (N-m-sec) 


h (N-m-sec) 


3 CMG/3 RW Combined Momentum Space 3 CMG/3 RW Combined Momentum Space 


3 te ■ ■ ft 

_ ■ - • ., 'f \ * . v j . • j? y. f,. ... . 




Figure 33. 


h x (N-m-sec) 

Composite array momentum volume for a 
RW/three CMG hybrid MED system 


three 
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B. TORQUE ENVELOPES 

1. Reaction Wheels 

With the momentum vector for the RW being essentially 
fixed along the rotation axis (assuming a perfectly 
balanced wheel), the torque space assumes the same shape as 
the momentum space, as seen in the torque space for the 
four RW pyramidal array in Figure 34. The other 

configurations are similar to the example shown in the 
previous section. 


Four RW Torque Space 



Figure 34. Four RW pyramid array maximum torque space 

2. CMG Torque Space Analysis 

a. Reduced Dimensional State Space R 3 Analysis 

When dealing with state space representations of 
a CMG array, as discussed in Chapter III Section C, it is 
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very important to maintain the contiguous or semi- 
contiguous condition when reducing dimensionality. 

There are several possible metrics that can be 
calculated for the state space analysis. One typical 
analysis is based on testing the singularity condition (SC) 
of a particular state as defined by Equation (102) [8, pg. 

467] . 

SC = det(A4 r ) (102) 

Where A can be defined as in Equation (103). 



For the case n> 3, a non-zero SC implies that the 

pseudoinverse control allocation, A + = A r ^AA r ) , exists 

because AA r is invertible. Thus, a vector of gimbal rates 
can be found that will satisfy Equation (103) . Therefore, 
developing a state space map based on the SC of each state 
allows a region of gimbal angles for each CMG to be created 
that is singularity free. 

Figure 35 shows a map of the SC based on state 
index value, with Figure 36 showing the gimbal angles based 
on state index number with zero representing the bottom 
left corner of the singularity free box and increasing to 
the right and up for state index one and two, respectively. 
The state indexes themselves were created by taking pairs 
of CMGs (one and two for the first state index, three and 
four for the second) and allowing the gimbal angle of one 
CMG in each pair to cycle through an entire circle before 
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incrementing the other CMG in the pair. This continued 
until the second CMG gimbal angle had gone through 360° . 
Figure 37 describes the process in pseudo code. 


Raw SC for 4 CMG Pyramid 
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Figure 35. 


Four CMG 
space 


pyramid reduced dimensional state 
analysis showing raw SC data 
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Allowable Gimbal Angles 



Figure 36. Four CMG pyramid reduced dimensional state 

analysis gimbal angle constraints 


forgimbal2 = 1 to 360 
forgimball = 1 to 360 
stateindexl(kk) = [gimball,gimbal2] 
gimball = gimball + increment 
kk = kk + 1 
end 

gimbal2 = gimbal2 + increment 
end 

Figure 37. Pseudo code representation of the creation 

of the state indices 


Figure 38 shows the state conditions that meet 
the inequality SC > 0. The black box present in Figure 36 
and Figure 38 indicate the identified singularity free 
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region, based on a simple search algorithm that finds the 
largest square area that meets the required condition. 


220 

200 

180 

160 

CM 

X 140 

m 

c 120 

£ ioo 

JS 

w 80 
60 
40 
20 


Figure 38. Four CMG pyramid reduced dimensional state 

analysis showing analyzed SC data 

Although the gimbal limits can be pulled off of 
Figure 37, understanding that the steps seen in CMG four 
represent the resolution of gimbal angles analyzed 
underscores the fact that this representation is not very 
accurate. The purpose of using such large steps initially 
is to reduce the analyzed state space since each gimbal has 
to be varied through the entire 360°. By using the 
resulting gimbal limits in a subsequent iteration and 
maintaining the number of steps for each angle range, the 
results can be refined at the same speed as the initial 
examination. This can be continued until the desired 
accuracy is developed. Figure 39 shows the results for the 
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second and third iteration while Table 2 gives the gimbal 
angle limits for each iteration. 


Iteration 


CMG 1 
(deg) 

CMG 2 
(deg) 

CMG 3 
(deg) 

CMG 4 
(deg) 

1 

Low 

0.00 

77.14 

0.00 

77.14 

High 

360.00 

180.00 

360.00 

180.00 

2 

Low 

0.00 

57.30 

0.00 

155.52 

High 

360.00 

118.68 

360.00 

229.18 

3 

Low 

0.00 

65.48 

0.00 

139.15 

High 

360.00 

122.78 

360.00 

196.44 


Table 2. Gimbal limits per iteration for four CMG pyramid 
reduced dimensional state space analysis based on 

SC 


The point-cloud momentum space resulting from the 
calculated gimbal angle limits is given in Figure 40. For 
these gimbal angles, there is a much larger proportion of 
volume in the positive x-y-z octant. So although the 
momentum space volume may be sufficient, the location of 
this volume with respect to the origin may not. For 
comparison purposes. Figure 41 shows the limited momentum 
space plotted within the unlimited momentum space. 
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Iteration 3 
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Figure 39. 


Results of iteration two and three for the four 
CMG pyramid reduced dimensional state space 
analysis based on SC 
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Gimbal Limited Momentum Space 



h x (N-m-sec) h x (N-m-sec) 

Figure 40. Momentum space for gimbal angle limited four CMG 
pyramid reduced dimensional state space analysis 

based on SC 

A useful, though not exhaustive, test to 
determine if the volume is singularity free is to pick two 
points within the volume and use a simple pseudoinverse 
controller to drive the four CMG array from one momentum 
state to another. Under normal circumstances, the 

pseudoinverse control law tends to drive the array towards 
a singularity, but by staying within the declared 
boundaries singularities can be avoided. 
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Figure 42 and Figure 43 show 
test, using a simple PD controller to 
error to zero. 


the result of this 
drive the momentum 




Gimbal Limited Momentum Space Comparison 


Gimbal Limited Momentum Space Comparison Gimbal Limited Momentum Space Comparison 


-sos 

h x (N-m-sec) 


-sos 

h x (N-m-sec) 


Figure 41. Reduced dimensional limited- and unlimited- 
gimbal angle four CMG pyramid momentum 
space comparison 

Normal attitude control operations do not 
calculate actual angular momentum values; they are usually 
based on current and final quaternions. If doing a 
standard eigenaxis maneuver using quaternions, the 
intermediate step of calculating angular momentum and 
ensuring the trajectory stays within the bounds (difficult 
due to the irregular shape of the volume) can be neglected 
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by simply ensuring all CMGs stay within their gimbal angle 
bounds. 


Momentum Error vs Time 




Torque vs Time 



Singularity Condition vs Time 



Figure 42. Pseudo-inverse control law operation in 
singularity free CMG momentum space 


When using the pseudo-inverse control law as a 
basis, there are several difficulties to overcome when 
designing this new implementation. First of all, 
controlling gimbal rate behavior near gimbal angle 
boundaries is very important. It would be very undesirable 
to generate a high gimbal rate in an adverse direction, 
forcing the CMG across the operational boundary. To 
minimize rate changes, one can simply convert to a delta 
form of the control law as given in Equation (104) . 
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( 104 ) 


8 = AS + 8 0 
h = A^AS + <5q) 

A S = A + (ii-AS 0 ) 


Gimbal Limited Momentum Space 




-5 -4 -3 -2 -1 o 1 2 3 

h x (N-m-sec) 
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Gimbal Limited Momentum Space 


Figure 43. Pseudo-inverse control 

singularity free 


h x (N-m-sec) 

Gimbal Limited Momentum Space 



h x (N-m-sec) 


law momentum trajectory 
CMG momentum space 


in 


The properties of the least squares solution will 
then minimize |Aj| , thus minimizing gimbal rate changes 

[16] . 

It is also necessary to ensure that the gimbal 
rates tend to avoid the saturation conditions. To make 
this happen, a weighting factor, w f , can be added to 
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Equation (94) . This factor will act to drive the solved A 8 
such that AS + S 0 will direct the gimbal away from the limit. 
The full steering control law is given in Equation (105) . 


h = AS + <5 0 + w f 'j 
A 8 = A + (h-AS + Aw f \ 


(105) 


Although no analysis has been completed to 
determine create a working example of Equation (105), it is 
likely that the design of a w f function to achieve the 
desired behavior will be the most challenging part. 


b. Full Dimensional State Space R' ,+1 Analysis 

One of the primary benefits of reduced 

dimensional analysis lies with the fact that intermediate 
results can be represented in a form that is fully 
representable in R 3 . A close inspection of the state index 
plot in Figure 39 shows that not all permutations of CMG 
gimbal angle states are analyzed. Though the states that 

are analyzed are contiguous, there are additional 
contiguous states that are ignored. This is the heart of 
why the previous analysis was characterized as semi- 
contiguous. The unanalyzed states impose unnecessarily 

restrictive constraints on the CMG array based on the 
gimbal angle limits. 

A fully contiguous analysis representation 

requires a five dimensional plot for the system being 
analyzed. This plot cannot be easily displayed for 
visualization with current technology. Fortunately, the 
same analytical principles still apply. 
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By storing the SC values in a four dimensional 
matrix, a hypercube can be formed within that matrix to 
determine the gimbal angle limits for singularity free 
space. With these limits, similar analyses such as the 
momentum space point cloud can be performed. 

Table 3 provides the results for each iteration 
of the full dimensional state space analysis while Figure 
44 and Figure 45 give the momentum space and provide a 
comparison to the prior methods, respectively. 


Iteration 


CMG 1 
(deg) 

CMG 2 
(deg) 

CMG 3 
(deg) 

CMG 4 
(deg) 

1 

Low 

25.71 

102.86 

25.71 

102.86 

High 

180.00 

257.14 

180.00 

257.14 

2 

Low 

0.00 

193.40 

0.00 

193.40 

High 

113.65 

322.48 

113.65 

322.48 

3 

Low 

0.00 

167.90 

21.51 

199.92 

High 

86.05 

295.97 

107.57 

327.98 


Table 3. Iterative results for the full dimensional state 

space analysis of the four CMG pyramid 

The complexities of the gimbal limited shapes 

make comparison between the two methods difficult. One 
fairly simple metric to compare the full and reduced 

dimensional shapes is to compare the radii of the largest 
sphere centered at the origin that can be inscribed within 
the volume of the shapes. For the reduced dimensional 
analysis case, this radius is approximately 0.1241 h N-m-s 
while for the full dimensional analysis case it is 

approximately 0.0689 h N-m-s. The full dimensional 

momentum space, however, is much more centered than the 
reduced dimensional analysis momentum space. The spheres 
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are shown superimposed on the respective shapes in Figure 
46 and Figure 47. 


Gimbal Limited Momentum Space 




Gimbal Limited Momentum Space 


Gimbal Limited Momentum Space 


-2-10123 

h x (N-m-sec) 



Figure 44. 


Four CMG pyramid full dimensional state space 
analysis-based singularity free momentum space 


Although Crews ([4]) indicates that the maximum 
singularity-free sphere has a radius of 1.0069 h N-m-s, a 
number which greatly exceeds the normalized radius of 
0.0689 h N-m-s for the full dimensional state space 

analysis and 0.1241 h N-m-s for the reduced dimensional 

state space analysis, the difference 
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is expected. 


The 

























1.0069 h N-m-s sphere is based only on neglecting elliptic 

singularities while this state space analysis avoids all 
singularities. Reanalyzing the gimbal limits by allowing 
hyperbolic singularities within the operating momentum 
volume will increase the available performance. 
Additionally, the gimbal limits are based on finding the 
largest singularity-free square or hypercube due to the 
simplicity of implementing the search algorithm as well as 
the simplicity of analyzing the resulting gimbal limits. 
This provides an additional performance constraint as it 
does not include all contiguous states in the maximum area 
or hypervolume that are singularity free. Finally, a voxel 
based maximal surface detection algorithm may improve the 
performance of the state space analysis with respect to the 
maximum sphere metric. 
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Gimbal Limited Momentum Space Comparison 




Gimbal Limited Momentum Space Comparison Gimbal Limited Momentum Space Comparison 
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Figure 45. Four CMG singularity free momentum space 

comparison 

Due to the difficulty in visually comparing the 
complex surfaces of singularities [4,17] and the complex 
shape of Figure 44, further study is required to verify 
that the results of the state space analysis matches those 
found in the literature. 

Additional work is required to fully determine 
whether reduced dimensional state space analysis provides 
sufficient performance or if a full dimensional state space 
analysis is required to maximize array performance. 
Ultimately, the result may depend on the application. 
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Gimbal Limited Momentum Space 
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Figure 46. Maximum sphere superimposed on the reduced 

dimensional state space analysis-based singularity 

free momentum space 
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Gimbal Limited Momentum Space 



h x (N-m-sec) 

Gimbal Limited Momentum Space 



h (N-m-sec) -2-10123 

h x (N-m-sec) 

Figure 47. Maximum sphere superimposed on the full 

dimensional state space analysis-based singularity 

free momentum space 


C. REACTION FORCES 

1. Validation and Verification 

Due to the lack of material in literature as compared 
to torque and momentum space studies, rigorous validation 
and verification of the reaction force equations was 
required. The general methodology used to conduct the 
validation and verification involved identifying particular 
initial conditions that output a simple, intuitive result, 
preferably involving a single parameter such as axial 
torque or transverse force. 
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There are two primary aspects that need to be 
verified. The first is the inertial motion of the MED. 
This is the right hand side of the force-balance equation. 
The inertial motion equation determines the linear and 
rotational trajectories based on known system geometry and 
simulated or real-world telemetry. The other aspect is the 
left hand side of the force-balance equation. It is the 
summation of all moments and forces and is used to 
determine the reaction forces and moments. Within those 
two primary sections are a total of six parts to be 
verified: forces in the axial and both transverse 

directions and moments in the axial and both transverse 
directions. 

a. Reaction Wheel Inertial Motion Verification 

For the reaction wheel tests, the initial 
conditions and results are given in Table 4 through Table 
9. The RW is assumed to be oriented such that the spin 
axis is parallel with the z-axis of the spacecraft and has 
the following physical characteristics: 

'1 0 0 " 

I m = 0 1 0 kg-nr 

0 0 1 

m = 1 kg 

The first test examines rotations in the axial 

direction and is broken up into four parts. The first two 

parts have the RW centered in the spacecraft body and have 

the RW accelerate at a constant rate and the spacecraft 

accelerate at a constant rate, respectively. The results 

for these tests give the expected inertial moment in the 

axial direction. The third part involves the RW 
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accelerating at a constant rate in one direction and the 
spacecraft accelerating at the same rate in the opposite 
direction. These two motions cancel each other out, 
resulting in no inertial moments or forces applied to the 
RW. The final part has the RW displaced from the center of 
the spacecraft and has the spacecraft rotating at a 
constant rate. The result shows an inertial force oriented 
towards the center of the spacecraft. 



Applied Ang. Accel. 
(rad/sec 2 ) 

Initial Ang. Rates (rad/sec) 


RW Location (m) 

S/C 

RW 

S/C 

RW 

Xs/c 

Ys/c 

^ 3 !! 

X 

y 

z 

z 

fl 

n 

z 

z 

Test 1 

a 

0 

0 

0 

0 

0 

0 

0.1 

0 

0 

0 

0 

b 

0 

0 

0 

0 

0 

0.1 

0 

0 

0 

0 

0 

c 

0 

0 

0 

0 

0 

- 0.1 

0.1 

0 

0 

0 

0 

d 

0.1 

0 

0 

0 

0 

0 

0 

0 

0 

a 

0 


Table 4. RW V&V test one initial conditions 



RW Inertial Force (N) 

RW Inertial Moment (N-m) 


X 

KM 

z 

X 

y 

z 


a 

0 

0 

0 

0 

0 

0.1 

1T 

0 

0 

0 

0 

0 

0.1 

c 

0 

0 

0 

0 

0 

0 

~d” 

- 0.1 

0 

0 

0 

0 

0 


Table 5. RW V&V test one results 



Applied Ang. Accel. 
(rad/sec 2 ) 

Initial Ang. Rates (rad/sec) 



RW COM Location (m) 

S/C 

RW 

S/C 

RW 






X 

y 

z 

z 

D 

n 

z 

z 

CM 

a 

0 

0 

0 

0 

0.1 

0 

0 

0 

0 

0 

0 

+■* 

<Sl 

0) 

b 

0 

0 

0 

0.1 

0 

0 

0 

0 

0 

0 

0 

H 

c 

0 

0 

0 

0 

0.1 

0 

0 

0 

0 

0 

0.2 


Table 6. RW V&V test two initial conditions 
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RW Inertial Force (N) 

RW Inertial Moment (N-m) 


X 

y 

z 

X 

y 

z 

CM 

a 

0 

0 

0 

0 

0.1 

0 

+■* 

CO 

CL) 

~V 

0 

0 

0 

0.1 

0 

0 

H 

c 

0 

0 

0 

* 

0.1 

0 


Table 7. RW V&V test two results 


* 


01 

T 


/ (iV - m ) 


The second test is composed of three parts and 

examines transverse motions. For all parts of this test, 
the RW is located at the center of the spacecraft. Parts 
one and two have the spacecraft with a rotational 
acceleration along a single transverse axis. The results 
show an inertial moment along the corresponding axis with 
the appropriate magnitude. The third part repeats the 
first, however the RW is now rotating at a constant rate 

along its spin axis. The results show a constant inertial 

moment along the appropriate transverse axis as well as a 
linearly increasing inertial moment about the other 

transverse axis. This increasing moment matches the 

expected behavior of a CMG, which is what the initial 
conditions for the third part simulated. 



Applied Linear Accel, (m/sec 2 ) 



RW COM Location (m) 

S/C 



Xs/c 

Ys/c 

Zs/c 

X 

y 

z 

m 

a 

0 

0 

0 

0.1 

0 

0 

+■» 

CO 

0) 

b 

0 

0 

0 

0 

0.1 

0 

H 

c 

0 

0 

0 

0 

0 

0.1 


Table 8. RW V&V test three initial conditions 
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RW Inertial Force (N) 

RW Inertial Moment (N-m) 


X 

y 

z 

X 

y 

z 

m 

a 

0.1 

0 

0 

0 

0 

0 

+■* 

CO 

CL) 

~V 

0 

0.1 

0 

0 

0 

0 

H 

c 

0 

0 

0.1 

0 

0 

0 


Table 9. RW V&V test three results 


The final test involved linear acceleration along 
all three RW axes. The results all showed an inertial 
force with the appropriate magnitude in the appropriate 
direction. 


b. RW Transverse Reaction Verification 

The reaction verification is broken up into two 
parts: analyzing an applied inertial transverse force and 
an applied inertial transverse moment. Both the moment and 
force are applied to act on the center of mass (COM). Each 
part is then subdivided into four sections. The first 
section analyzes the effect of moving the COM from one 
bearing to the other, along the length of the shaft. The 
second section analyzes the effect of shaft length. Both 
sections are performed for both the cantilever and fixed- 
fixed end condition cases. The RW is assumed to have the 
same physical properties as described in the inertial 
motion verification. Additionally, the nominal shaft 
length is one meter and the nominal COM location is at the 
midpoint of the shaft. The applied force has a magnitude 
of 1 N and the applied moment has a magnitude of 1 N-m. 
These values were chosen to generate plots whose 
relationships could be scaled to fit the scale of the MED 
array being analyzed. 
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Figure 49 through Figure 52 give the transverse 
force case with Figure 48 giving the direction of the 
applied force. In accordance with the convention described 
in Chapter V Section B, bearing one is at the bottom of the 
RW and bearing two is at the top. 
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Figure 48. Applied inertial force geometry for 

RW reaction validation 
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COM Location Variation 




Figure 49. Fixed-fixed RW COM variation reaction 
validation for an applied inertial force 
at the center of the shaft 
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Figure 50. Fixed-fixed RW shaft length variation 

reaction validation for an applied inertial 
force at the center of the shaft 
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COM Location Variation 
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Figure 51. Cantilever RW COM variation reaction 
validation for an applied inertial force 
at the center of the shaft 
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Figure 52. Cantilever RW shaft length variation 
reaction validation for an applied 
inertial force at the center of the shaft 
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For both the fixed-fixed and cantilever cases, 
the reaction forces and moments behave as expected with 
respect to the geometrical variations. 

Figure 53 shows the geometry for the applied 
moment section, with Figure 54 to Figure 57 giving the 
results. 


Applied Moment Geometry 
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Figure 53. Applied inertial moment geometry for 

RW reaction validation 
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COM Location Variation 




Figure 54. Fixed-fixed RW COM variation reaction 
validation for an applied inertial moment 
at the midpoint of the shaft 
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Figure 55. Fixed-fixed RW shaft length variation 

reaction validation for an applied inertial 
moment at the midpoint of the shaft 
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COM Location Variation 
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Figure 56. Cantilever RW COM variation reaction 
validation for an applied inertial moment 
at the midpoint of the shaft 
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Figure 57. Cantilever RW shaft length variation 

reaction validation for an applied inertial 
moment at the midpoint of the shaft 
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All reaction forces and moments acted as expected 
through this analysis, which validate the application of 
the beam tables. 

c. Reaction Wheel Axial Reaction Verification 

Due to the assumed construction of the fixed- 
fixed system, bearing one supports all axial force and 
moment loads. This matches the reaction of the cantilever 
system. Since Equations (71) and (81) explicitly enforce 
this property, additional analysis is not required. 

d. CMG Inertial Motion Verification 

As the inertial motion of the gimbal frame 

follows the same kinematics as the RW frame, the prior 
analysis sufficiently proves the equations used. The rotor 
frame, however, has the added complexity of gimbal motion 
and so must be analyzed. 

Due to the complexity of the resulting inertial 

forces and moments on the rotor from combined gimbal and 
spacecraft rotations, the utility of performing as 
exhaustive an analysis as was performed on the RW equations 
is limited. Therefore, two specific rotor cases will be 

examined: motion by the gimbal only and motion by the 
gimbal with counter-motion by the spacecraft (similar to 
test l.c for the RW) . The mechanical properties of the 

rotor will match those of the RW, the origin of the gimbal 
frame will be coincident with the origin of the spacecraft 
body frame, and the gimbal rotation axis will be parallel 
to the Z-axis of the spacecraft body. All bodies will 
start at rest. Table 10 gives the specific conditions of 
the test. 
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Applied Ang. Accel, (rad/sec 2 ) 


Rotor COM Location (m) 

S/C 

Gimbal 

Xg 

Y g 

Z G 

X 

y 

z 

8 

Test 1 

0 

0 

0 

0 

0 

0 

0.1 

Test 2 

0 

0 

0 

0 

0 

-0.1 

0.1 


Table 10. CMG V&V initial conditions 


Figure 58 and Figure 59 give the results of the 
tests. The small constant positive inertial force in the y 
direction in test one is caused by the angular acceleration 
of the gimbal. The decreasing force in the negative x 
direction represents the centrifugal force caused by the 
actual rotation of the gimbal. As the gimbal rate 
increases, the magnitude of the force increases as seen in 
Equation (106) . 

"«, = V, + ^([^] + [^-][^]) ;R g r r 

+2 ^[^ 6W ]^[®* /4 ]*r r + ^//?([g* /i ] + [®* /6 ][®* /6 ])*r r (10 6) 

N a r = N b R b g R([q g,b ] + [^ ,b ][a g,b ]) g r r 
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Inertial Motion CMG Test 1 
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Results of the first CMG V&V test 
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Figure 59. Results of the second CMG V&V test 
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For the second test, the prominent feature is the 
noise in the inertial force plot. Though of a negligible 
magnitude, it likely appears due to several adverse factors 
including the rudimentary simulation method and relatively 
rapidly rotating bodies. The important feature is the fact 
that the counter-spinning bodies remove the inertial moment 
from the rotor as it now appears to stand still in the 
inertial frame. 

2. Reaction Wheel Maneuvers 

The spacecraft geometry for the RW analysis is given 
in Figure 60 and tabulated in Table 11. The maneuver is an 
arbitrary rest to rest slew from one orientation to 
another. The analyzed telemetry data includes only the 
spacecraft quaternions and body rates, the RW angular 
acceleration and an array of the time steps. 


RW Array Layout 



— RW 1 
RW 2 

— RW 3 

— RW 4 

* BRG1 
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* COM 
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Figure 60. Geometry detail for the RW maneuver force 

and moment analysis 
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RW 

r ( de §) 

P (deg) 

COM (X b /Y b /Zb) (m) 

1 

0 

54.73 

0.25 

0 

0 

2 

90 

54.73 

0 

0.25 

0 

3 

180 

54.73 

-0.25 

0 

0 

4 

270 

54.73 

0 

-0.25 

0 


Table 11. Tabulated values for RW array layout for 

maneuver analysis 

The RWs used in the analysis have the same properties 
as those used in the momentum space case study. Telemetry 
was derived from the model developed in [4]. 

a. Eigenaxis Maneuver 

The first maneuver is based on an eigenaxis 
controller where the eigenaxis between the initial and 
final state is calculated and a torque command is developed 
to rotate the spacecraft about this axis. Figure 61 
describes the basic telemetry for the maneuver. 


ill 




Wheel Rates 



Body Rate Norm 



Figure 61. Eigenaxis slew maneuver for RW array force 

and moment analysis 

During the maneuver, the first and third RW 
saturate at +/- 4000 RPM to achieve a constant spacecraft 
body rate magnitude of slightly above 4 deg/sec about the 
eigenaxis. The calculated inertial force and moment 
magnitudes are given in Figure 62 and Figure 63. 
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Inertial Force Magnitude 



Figure 62. Inertial force magnitudes during an eigenaxis 

slew maneuver 


Inertial Moment Magnitude 



time (sec) 


Figure 63. Inertial moment magnitudes during an eigenaxis 

slew maneuver 


During this maneuver, the magnitudes of the 
inertial forces are equal for opposite pairs of RWs (one 
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and three, two and four). This matches expectations as the 
RWs lie equidistant along a common axis from the spacecraft 
body origin which is the point about which the spacecraft 
rotates. Thus, one would expect the centrifugal force on 
the RWs to be equal in magnitude and opposite in direction. 
The difference in magnitude between the two pairs indicates 
that RWs two and four are moving faster through inertial 
space, indicating that the projection of the eigenaxis 
along the line connecting RW two and RW four is larger than 
the projection on the line connecting the other set. This 
assumption is borne out by the magnitude of the inertial 
moments. Since the spacecraft is spinning about an axis 
that is closer to the rotation axis of RW one and three, 
they will have a larger inertial moment induced on them. 

The steps seen at about 12 seconds and 46 seconds 
indicate the times when the RW angular acceleration goes 
from non-zero to zero (and vice-versa) . 

The reaction force plots for the fixed-fixed 
configuration given in Figure 64 and Figure 65 show an 
interesting property. The inertial moments are between one 
and two orders of magnitude larger than the forces. From 
Figure 54, a one meter shaft results in reaction forces 
about 1.5 times larger than the applied moment. 
Furthermore, Figure 55 shows that a 0.2 m shaft, the 
assumption used for this analysis, results in reaction 
forces eight times larger than the applied moment. The 
inertial moment, therefore, contributes most to the 
reaction force at the bearings. 
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BRG1 Rxn Forces 



Figure 64. Bearing one reaction force magnitudes for the 
eigenaxis slew, fixed-fixed configuration 


BRG2 Rxn Forces 



Figure 65. Bearing two reaction force magnitudes for the 
eigenaxis slew, fixed-fixed configuration 


That fact explains why both bearings have 
approximately the same force magnitude, despite the fact 
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that bearing two does not support a force in the axial 
direction. When examining the reaction moments. Figure 66 
and Figure 67, bearing one shows a significant larger 
moment than bearing two during wheel operation. This is 
attributed to the fact that only bearing one resists an 
axial moment due to the motor being co-located with bearing 
one. The larger reaction moment magnitude for RWs one and 
three stems from the fact that they are accelerating at a 
higher rate. 

When the RWs are not accelerat ing, the reaction 
moment magnitudes on both bearings roughly match each 
other. 


BRG1 Rxn Moments 



Figure 66. Bearing one reaction moment magnitudes for the 
eigenaxis slew, fixed-fixed configuration 
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BRG2 Rxn Moments 



Figure 67. Bearing two reaction moment magnitudes for the 
eigenaxis slew, fixed-fixed configuration 

With the cantilever case, the bearing one 
reaction forces and moments act as expected, shown in 
Figure 68 and Figure 69. The small inertial force 
magnitude and short lever arm prevent the creation of a 
noticeable reaction moment at the bearing, thus the 
reaction moment and inertial moment plots appear to match. 

The magnitude of the reaction force on bearing 
one in the cantilever case is much less than for the fixed- 
fixed case because, according to the standard beam tables 
[13, pp. 208-223], the inertial moment does not result in a 
reaction force on the bearings. 
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BRG1 Rxn Forces 



Figure 68. Bearing one reaction force magnitudes for the 
eigenaxis slew, cantilever configuration 
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Figure 69. Bearing one reaction moment magnitudes for the 
eigenaxis slew, cantilever configuration 


118 






























































b. Optimal Maneuver 


The second maneuver analyzed is an optimal 
maneuver where the cost function, J , represents a quantity 
that the controller minimizes. This cost function can 
include such quantities as elapsed time or wheel power. 
The same RW configuration is used and the telemetry is 
derived from the same source. The maneuver is conducted 
with the same initial and final conditions as before. 
Figure 70 shows the maneuver. 


Wheel Rates 



Figure 70. Optimal slew maneuver for RW array force and 

moment analysis 

For this maneuver, the wheel rates of all wheels are 
now closer to or at saturation for a longer period of time. 
The body rate magnitude is somewhat higher throughout the 
maneuver, and the final state is reached slightly sooner 
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than the eigenaxis maneuver. This additional exercising of 
the RWs is typical of an optimal maneuver as the algorithm 
seeks to use the excess performance that is present but 
that the eigenaxis slew ignores. 

Since there is no single axis of rotation for the 
maneuver and the RWs tend to continuously accelerate or 
decelerate, the inertial force and moment magnitudes 
induced on the RWs tend to be more complicated, as shown in 
Figure 71 and Figure 72. 


Inertial Force Magnitude 



Figure 71. Inertial force magnitudes during an optimal 

slew maneuver 


120 



























Inertial Moment Magnitude 



Figure 72. Inertial moment magnitudes during an optimal 

slew maneuver 

There are some similarities between the time optimal 
and eigenaxis slew maneuvers with regards to the inertial 
forces and moments. As with the eigenaxis slew, the 
optimal slew produces inertial forces that are equal in 
magnitude to the same paired RWs: one and three, two and 
four. This equality is expected due to the positional 
symmetry in the system. Since all RWs are equidistant from 
the z-axis, a spacecraft body rotation about this axis 
produces equal inertial force in all components. RW one 
and three are positioned along the x-axis equidistant from 
the origin, so a spacecraft body rotation rate about the y- 
axis produces equal and opposite inertial forces in that 
pair, and zero inertial force in the other. The opposite 
is true for a spacecraft body rotation rate about the x- 
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axis. The inertial forces tend to be slightly larger for 
all four components due to the higher body rates; however 
they are still an order of magnitude below the inertial 
moments. 

When analyzing the behavior of the moments, though, it 
is possible to determine when changes in wheel acceleration 
occur as these accelerations cause moments that tend to 
dominate. This is expected as the wheel can change 
acceleration much faster than the spacecraft will alter its 
rotational motion. Removing the torque induced by the 
wheel itself would produce much smoother plots, as seen 
during the periods of constant wheel acceleration in Figure 
63. 

The bearing force and moment reaction magnitudes are 
given in Figure 73 through Figure 76. As before, the 
inertial moments dominate the reaction forces, and bearing 
one supports much larger reaction moments than bearing two 
due to the actual torque produced by wheel accelerations. 
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BRG1 Rxn Forces 



Figure 73. Bearing one reaction force magnitudes for an 
optimal slew, fixed-fixed configuration 


BRG2 Rxn Forces 



Figure 74. Bearing two reaction force magnitudes for an 
optimal slew, fixed-fixed configuration 
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BRG1 Rxn Moments 



5. Bearing one reaction moment magnitudes for an 
optimal slew, fixed-fixed configuration 


BRG2 Rxn Moments 



Figure 76. Bearing two reaction moment magnitudes 

for an optimal slew, fixed-fixed configuration 
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Figure 77 and Figure 78 give the cantilever results. 
The reaction forces and moments match the inertial forces 
and moments, an expected result. 

Although the optimal maneuver results in more active 
RW utilization, the force and moment magnitudes seen by the 
bearings are of roughly equivalent magnitudes to those seen 
during the eigenaxis maneuver. Because the wheel 
accelerations tend to dominate the inertial moment equation 
(and the inertial moment dominates the inertial force with 
respect to bearing reactions), the optimal maneuver tends 
to maintain a higher average loading. This particular 
eigenaxis maneuver, however, produces very large and rapid 
changes in force and moment magnitudes due to the bang-off- 
bang nature of the maneuver. This shock to the system is 
not present in the optimal maneuver. 


BRG1 Rxn Forces 



Figure 77. Bearing one reaction force magnitudes for an 
optimal slew, cantilever configuration 
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BRG1 Rxn Moments 



Figure 78. Bearing one reaction moment magnitudes for an 

optimal slew, cantilever configuration 

3. CMG Maneuvers 

Using four CMGs arranged in the same geometrical shape 
shown in Figure 60 and described in Table 11, reaction 
analyses for an eigenaxis slew and an optimal slew were 
performed. In this case, the initial and final conditions 
between the two maneuvers are slightly different because 
the telemetry data is from a higher fidelity simulation of 
an imaging spacecraft. The CMG mechanical properties are 
the same as those for the momentum and torque space 
analysis in Chapter VI Sections A and B with the addition 
that the rotor shaft length is assumed to be 0.1 m, the 
gimbal shaft length is assumed to be 0.5 m and bearing one 
for the rotor is assumed to be mounted 0.05 m from the 
gimbal COM in the rotor's axial direction {x =x r ). The CMG 
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has an assumed mass of 2 kg. The assumed rotor inertia 
tensor is: 

".0429 0 0 

I RW = 0 .0429 0 kg ■ m 1 . 

0 0 .0025 

Although there are four possible variations for the 
gimbal and rotor end point conditions, only the fixed- 
fixed/fixed-fixed case will be examined as it provides 
sufficient detail for this analysis. A more rigorous 

examination would be warranted in other cases such as 
actual satellite or CMG mechanism design. 

a. Eigenaxis Maneuver 

The analyzed maneuver is actually the end of one 
maneuver and the beginning of another. Due to the nature 
of the data, representing an entire maneuver is impractical 
due to the length of time involved. Analyzing the impulses 
between two maneuvers reduces the amount of time to be 
analyzed while still examining two major control events. 

Initially the spacecraft is rotating at a small 
rate. The CMGs are gimbaled to reduce the body rate for a 
few seconds. They are then gimbaled again to induce a new 
spacecraft body rate. Figure 79 through Figure 85 give the 
maneuver telemetry and results of the force analysis. The 
large force of 36 N seen at 41 and 45 seconds in the first 
gimbal bearing is caused by the inertial moment induced on 
the rotor at those times. Due to the complexity of the 
system dynamics, this result cannot be easily gleaned from 
the analysis of the equations of motion as was done in 
Equation (106) . 
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Due to the short length of the rotor shaft, the 
reaction forces are many times larger than the inertial 
moment would imply. The reaction force at the rotor 
bearing one is then transmitted to the gimbal in the axial 
direction. Since the second gimbal bearing does not 

support forces in this direction, the entire load is taken 
up by the first bearing (or more specifically, the 
gimbaling motor located at that bearing). 



Body Rate Norm 



Figure 79. Eigenaxis slew telemetry for CMG force analysis 
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Figure 80. Inertial forces from the CMG force analysis 

eigenaxis slew 


Inertial Moment Magnitude 




Figure 81. Inertial moments from the CMG force analysis 

eigenaxis slew 
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BRG1 Rxn Forces 
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Figure 82. 


Bearing one reaction force magnitudes from the 

eigenaxis slew 


BRG2 Rxn Forces 




Figure 83. Bearing two reaction force magnitudes from the 

eigenaxis slew 


130 













































































BRG1 Rxn Moments 
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Figure 84. Bearing one reaction moment magnitudes from the 

eigenaxis slew 


BRG2 Rxn Moments 




Figure 85. 


Bearing two reaction moment magnitudes from the 

eigenaxis slew 
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Although the maneuvers between the RW and CMG 
analyses are very different, it appears that a CMG has the 
potential to induce much larger forces on its bearings 
compared to a RW. This is also how the CMG provides large 
torques to the spacecraft. Thus, the bearings for a CMG 
must be built to withstand larger loads than those for a 
comparably sized RW. 

b. Optimal Maneuver 

Similar to the eigenaxis maneuver, the optimal 
maneuver shows the completion and commencement of two back- 
to-back slews. The gimbal accelerations are not as 
excessive in this case which slightly reduces the induced 
rotor inertial moment. This, in turn, reduces the large 
reaction force seen at the gimbal/rotor junction and the 
first gimbal bearing. The results are given in Figure 86 
through Figure 92. 
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Figure 86. Optimal slew telemetry for CMG force analysis 
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Figure 87. Inertial force magnitudes from the CMG force 

analysis optimal slew 
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10 -3 Inertial Moment Magnitude 




Figure 88. Inertial moment magnitudes from the CMG force 

analysis optimal slew 
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Figure 89. Bearing one reaction force magnitudes for the 

optimal slew 
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Figure 90. Bearing two reaction force magnitudes for the 

optimal slew 
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Figure 91. Bearing one reaction moment magnitudes for the 

optimal slew 
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BRG2 Rxn Moments 




Figure 92. Bearing two reaction moment magnitudes for the 

optimal slew 


Contrary to what was seen for the RWs where the 
peak measurements between the control cases were roughly 
equal, the peak reaction forces for the eigenaxis maneuvers 
are slightly larger than those seen during the optimal 
maneuvers. The optimal maneuver, however, maintains the 
reaction force for a larger amount of time while the 
eigenaxis maneuver introduces more abrupt force changes. 
The same trend is seen for the reaction moments. This 
indicates that, for this particular maneuver with this 
particular geometry, the optimal controller is easier on 
the components than the eigenaxis maneuver. 
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D. SUMMARY 

Chapter VI presented a series of applications of the 
equations and algorithms derived in the prior chapters. 
Overall, the algorithms were shown to be able to rapidly 
analyze many different variations on a potential ACS, as 
well as the ability to analyze maneuvers to determine 
expected reaction forces and moments using only MED 
geometry and standard telemetry data, allowing for 
comparison of different controllers and different array 
geometries. 

1. Momentum Envelope 

For the RW momentum envelopes, the equations were 
validated by analyzing the geometry of the momentum spaces 
with the geometry of different arrays. Additionally, an 
analysis was performed on the effect of rotor imbalance and 
indicated that, for expected imbalance magnitudes, there 
was essentially no performance loss in the array. 

The CMG MIA analysis was effective in rapidly 
generating a highly detailed maximum momentum space for 
different CMG array geometries. The associated angle 
residual plots gave an indication where interior 
singularity surfaces met the maximum surface. 

Finally, the algorithms were used to generate a 
maximum momentum envelope for a hybrid RW/CMG array by 
summing the individually generated RW array momentum 
envelopes and CMG array envelopes. This resulted in the 
possibility of estimating the hybrid array performance with 
respect to spacecraft slews as well as introducing the idea 
of adding RWs to a CMG array to allow the CMGs to steer 
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allowing expanded 


around elliptic singularities thus 
utilization of the CMG performance envelope. 

2. Torque Envelopes 

The torque envelope was shown to be of the same shape 
as the momentum envelope for RW arrays. 

For CMG arrays, the state space analysis algorithms 
were implemented and used to generate singularity free 
momentum spaces. It was discovered that the current state 
space analysis methods resulted in a momentum envelope that 
allowed about 12% of the capability of current CMG array 
capability for the semi-adjacent analysis, and about 6% for 
the fully adjacent analysis. This had several causes, the 
primary one being that the current performance metric is 
based on including hyperbolic singularities while the state 
space analyses excluded all singularities. A theoretical 
control algorithm was presented, and an initial test was 
performed showing that both the control algorithm concept 
was valid and the momentum space held no singularities 
along the trajectory. 

3. Reaction Force Analysis 

The developed equations were validated for RWs and 
CMGs using simple test cases. Third party simulators were 
used to generate telemetry to analyze different 
controllers, analyze different RW and CMG construction 
assumptions and prove the feasibility of using third-party 
telemetry for the analysis. 

For both RW and CMG arrays, an eigenaxis slew and an 
optimal slew were analyzed. For the RW analyses, peak 

force and moment magnitudes were roughly equal. The 

138 



optimal slew operated the RWs at a higher duty cycle, but 
provided a much smoother reaction force and moment history. 
The eigenaxis slew tended to create abrupt changes in the 
reaction forces seen on the system. 

For the CMG arrays, the peak magnitudes for the 
optimal maneuver were less than for the eigenaxis maneuver. 
The smoothness seen in the RW analysis carried over to the 
CMGs for the optimal maneuver, and the abrupt applications 
of relatively large forces and moments was even larger for 
the CMG case than the RW case. 
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VII. CONCLUSIONS AND FUTURE WORK 


A. CONCLUSIONS 

The goal of this thesis was to develop the equations 
and algorithms necessary to automate the analysis of MED 
arrays for the purpose of ACS design and analysis. This 
involved looking at three specific areas: momentum envelope 
analysis, torque envelope analysis, and reaction force and 
moment analysis. 

For RWs, a simple matrix equation. Equation (23), was 
developed that can rapidly create a point cloud of momentum 
states. This provided a means to estimate possible ACS 
performance based on the minimum momentum point on the 
maximum surface. This, in turn, allowed array geometry or 
the physical characteristics of the RW (such as inertia 
tensor or maximum wheel speed) to be varied and the result 
of this variance with respect to array performance to be 
analyzed. It was also shown that the RW torque space has 

the same shape as the RW momentum space, with the torque 

space magnitude being defined by the maximum wheel 
accelerations. 

For CMGs, the MIA process was developed to rapidly 
analyze the maximum momentum surface by removing the need 
to perform a sweep analysis through all gimbal angle 
combinations. It also provided an indication of where 

interior singular surfaces meet the maximal momentum 

surface via an angle residual plot. 

When analyzing torque space, it was shown that CMG 

arrays produce torque that is state dependent. This led to 

the development of the state space analysis in order to 
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generate singularity free momentum spaces by limiting the 
gimbal angles of each CMG in the array. The idea of 
contiguousness was introduced to ensure that an enclosed 
area within the reduced dimensional state space analysis 
map would only contain states that yield a realizable 
gimbal angle trajectory. Both of the two primary state 
space analyses, reduced dimensional and full dimensional, 
successfully created singularity free spaces. 

The reaction force and moment analysis involved a 
topic that had received little discussion in the researched 
literature. The equations were developed using a 
combination of first principle derivations and standard 
reference tables. The analysis was formulated such that 
common telemetry data and known MED array geometry, 
consistent with the previous momentum and torque analyses, 
would be sufficient to solve the equations. Results from 
the verification and validation tests showed expected 
trends and provided useful relationships between MED 
geometry changes and resulting reaction force/moment 
changes. Initial telemetry testing indicated that optimal 
controllers and eigenaxis controllers created reaction 
forces and moments that were of the same magnitude for RWs, 
while for CMGs the optimal controller showed smaller force 
and moment magnitudes. Additionally, eigenaxis controllers 
tended to impact the system with sudden, large forces and 
moments while the optimal controllers resulted in smoother 
reaction forces and moments. In all cases, the average 
duty cycle of all MEDs operating with the optimal 
controllers was higher than the average duty cycle under 
the eigenaxis controllers. 
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Overall, this thesis showed that automated analyses 
for generic MED arrays was feasible, and should prove 
useful for the rapid iteration of initial ACS design. 


B. FUTURE WORK 


Although the goal of developing algorithms to rapidly 
analyze arrays of RWs and CMGs was achieved, there are 
several additional features that can be added to improve 
analysis fidelity and usefulness. 


1. Momentum Analysis 

• Very tight tolerances in fly-wheel balancing are 
required to ensure the long life of RWs and CMGs. 
This does allow for some imbalances, however, and 
could impact the performance of platforms that 
require very tight pointing requirements. 
Especially for RWs, it would be useful to modify 
Equation (28) so that it included tracking of RW 
frame rotation as the wheel rotates. 

• For CMGs, the rotor provides the bulk of the 
component's momentum. The rotation of the gimbal, 
however, adds some momentum and this should be 
included as part of the analysis. This may also 
benefit the singularity free state space analysis 
as it would then provide for the possibility of 
angular momentum vectors outside the plane of 
rotation. 

• The maximal surface detection algorithm is a very 
rudimentary one that sections a unit sphere and 
finds the point furthest from the point in any 
given section. Anomalies can present themselves 
at the north and south poles of the sphere. A 
voxel-based approach could provide a more robust 
solution as well as a possible performance 
improvement. 

• Though the momentum algorithm can easily be used 
to find the extreme vertices of a RW momentum 
envelope, identification of connected points must 
be done manually. Automating this task would 
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provide an extremely rapid creation of the edges 
of the geometry. 

• The generation of the quasi-maximum space for a 

hybrid array currently depends on a very 
simplistic method examining every CMG point, 
finding the nearest RW point and adding the 
momentum values, subject to an angle-error 
constraint. Due to the highly irregular 

distribution of the RW points and the very 
regular distribution of CMG points, there can 
often be a miss where a CMG point has no 
corresponding RW point. Due to the geometrically 
simple shape of RW momentum spaces, it should be 
possible to determine the momentum value at any 
point on the surface and add that value to the 
CMG point. 

• Further analysis should be conducted to determine 
the reasonableness of including small RWs with a 
CMG array to allow torque-free elliptic 
singularity avoidance within the entire CMG 
momentum space. 

2. Torque Analysis 

• The box/hypercube singularity free search 
algorithm misses a significant amount of adjacent 
singularity free space due to the simplicity of 
implementation for both the search algorithm and 
the resulting singularity free momentum space 
analysis. This unnecessarily restricts CMG 
capability. An improved search algorithm would 
restore some of the lost performance. 

• Due to the numerical nature of the singularity 
free space created by the state space analysis, 
another method must be developed to verify that 
the space is actually singularity free. Ideally, 
this could be done analytically as it would 
possibly allow an analytic determination of 
allowable gimbal angle ranges, obviating the need 
for a numerical analysis. 

• Adding limited gimbal angles to the momentum 
inversion algorithm would allow the singularity- 
free momentum space to be created with the same 
method as the normal maximal surface. 


144 



• Further development of the control law in 

Equation (94) would allow full utilization of the 
singularity free space provided by the gimbal 
angle limitations. 

3. Reaction Analysis 

• The value of the variables and parameters 

involved in the force analysis can vary across 
many orders of magnitude. To prevent the 
insertion of numerical errors when performing 

arithmetic between very large and very small 

numbers, scaling should be introduced to the 
equations. 

• The simplified model may not accurately model the 

physical MEDs. Further analysis should be 

performed to improve the simplifications. 

• The design of a MED array could be optimized for 

a specific mission by assuming the operating 
regime of a spacecraft and designing the MED 

based on the calculated loads 
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